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Abstract 

Wc propose an approach that views \J(N C ) Yang-Mills theory as the critical point of an induced 
gauge model on the lattice. Similar recent proposals based on the color-flavor transformation rely 
on taking the limit of an infinite number of infinitely heavy particles. In contrast, we couple a finite 
number N-^ of auxiliary boson flavors to the gauge field and argue that Yang-Mills theory is induced 
when iVb exceeds N c and the boson mass is lowered to a critical point. Using the notion of Howe 
duality wc transform the induced gauge model to a dual formulation in terms of local gauge invariant 
variables. In the abelian case the Howe duality transform turns out to coincide with the standard one, 
taking weakly coupled U(l)<2=4 to strongly coupled 7Ld=^ lattice gauge theory. 

1 Introduction 

In what has come to be called "induced QCD", one starts from a theory of auxiliary fields 
coupled minimally to a gauge field background with gauge group G, the notable feature being 
that the gauge field by itself has no dynamics. A dynamical theory of gauge fields, typically 
with nonlocal gauge interactions, is induced by elimination of the auxiliary fields. Such 
an approach El El E] has been suggested on the lattice and in the continuum, using scalar 
fields as well as fermions in the fundamental or adjoint representations of G. Perhaps the best 
known model in the induced QCD category is the Kazakov-Migdal model [I], which attracted 
a flurry of interest in 1992-93, presumably because it admits a solution in the large-iV limit. 
There exist two recent papers El motivated by the color-flavor transformation that revive 
the old idea of inducing QCD on the lattice. 

The big question looming over all these induced gauge models is whether they do indeed lie 
in the universality class of Yang-Mills theory with gauge group G as intended. In most of the 
past work, the concern was with nonlocality of the effective gauge field action. To suppress 
nonlocal contributions, one used the trick of sending both the number of auxiliary fields and 
their mass to infinity in a specific manner. In the present paper, we propose a rather different 
approach. We choose the dynamics of the auxiliary fields to be local to begin with, so that 
the induced gauge field action automatically has that property. Then, fixing the number N\> 
of auxiliary fields, we tune the mass to a critical point so as to induce Yang-Mills theory. In 
order for this to work, it is crucial that bosonic scalars (as opposed to fermions) be used. 
In more detail, the effective action of the induced gauge model we propose is a sum over 
elementary plaquettes on a d-dimensional lattice: S = —2N^ ^ p Re Tr In (m — U(dp)) , with 



mass parameter m > 1. At unit mass the model has a critical point with a diverging corre- 
lation length, allowing a continuum limit to be taken. We claim that the m — > 1 continuum 
model flows under renormalization to Yang-Mills theory, provided that exceeds a certain 
minimum value. (What that minimum value is depends on the type of gauge group, its rank, 
and the space-time dimension.) The mathematical basis for this claim is that the local weight 
function U i— > |Det(m — U)\~ 2Nh on a classical compact Lie group G approaches for m — > 1 
the 5-function supported at U = 1, if N\, is large enough. 

The second major theme of the present paper is a duality transformation for the induced 
gauge model. While "duality" has been very successful for gauge theories in the continuum 
and with supersymmetries — important examples are electric-magnetic (or Montonen-Olive) 
duality, which constitutes an important aspect of the Seiberg-Witten solution for the low- 
energy dynamics of J\f = 2 super symmetric gauge theory with gauge group SU(2); or the 
Maldacena conjecture (or AdS/CFT-duality) , by which N = 4 super Yang-Mills theory 
with gauge group SU(iV) is believed to be dual to a type II-B superstring theory — , our 
interest here will be in pure gauge theories on the lattice, and without supersymmetry. In 
this class of quantum field theories duality is a well developed concept only for the abelian 
case [3 El El EH EEH EG] , where passing to the dual theory essentially amounts to taking a 
Fourier transform. 

A direct transcription of the abelian duality transform to nonabelian lattice gauge theories 
makes use of a character expansion of the plaquette statistical weight function (see f° r 
a recent reference). Integration over the gauge field then produces sums over products of 
Racah coefficients of G [these are higher-rank generalizations of tensor invariants known as 
3j, 6j, 9j, 12j symbols etc. for the case of G = SU(2)]. The semiclassical asymptotics of such 
coefficients is not easy to handle, and hence the continuum limit of the theory in this dual 
formulation remains unclear except in some special situations. 

We consider it to be an interesting feature of the determinant-type models we are going 
to introduce, that an alternative approach is possible. Our induced gauge models can be 
transformed to a dual description — for any one of the classical compact gauge groups G = 
XJ(N), Sp(2iV), and O(N) — by viewing the determinant weight functions as traces in an 
(auxiliary) Fock space. The gauge group G acts in this Fock space as one member of a dual 
pair of Lie groups, where "dual" is meant in the sense of R. Howe |14| 115]. Integration over 
the gauge fields simply projects on the G invariant subspace of Fock space. The other member 
of the Howe dual pair is a noncompact Lie group acting irreducibly on that subspace. This 
allows a fairly transparent description of what the dual lattice theory is (although we do not 
yet understand its continuum limit). For G = U(l) our Howe duality transform reproduces 
standard abelian duality upon elimination of some redundant degrees of freedom. 
Howe pairs underlie the color-flavor transformation |16] ITT] , which was originally conceived 
in the context of random matrix theory and disordered electron systems, and has recently 
been applied to the strong coupling limit of QCD ^1 ^3 HOI US • Howe pairs have also been 
used recently for the investigation of determinant correlations for quantum maps (221 • 

This paper is organized as follows. In Section |2] we place (local) bosons and/or fermions on 
the sites of an arbitary lattice and couple them to gauge fields in the standard manner to 
induce a pure lattice gauge theory with gauge group G. For the case of G = U(iV c ), Section^] 
establishes the <5-function property of the induced plaquette weight function when the number 
of boson flavors N\> > N c . The resulting boson induced gauge model is subjected to a careful 
investigation in d = 1 + 1 dimensions in Section We show that its partition function in the 
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continuum limit (and on any Riemann surface) agrees with the known partition function of 
some U(iV c ) Yang-Mills theory if > N c + 1. For = N c we identify an exotic continuum 
gauge theory, which is not Yang-Mills but of Cauchy type. Then, after introducing the notion 
of Howe pairs in the two-dimensional context, we subject the boson induced gauge model to 
a duality transformation in any space-time dimension in Section El We conclude with a 
summary and listing of open questions in Section H3 commenting in particular on the issue of 
convergence of the color-flavor transformation when the number of boson flavors is large. 

2 Fermion and Boson Induced Gauge Model 

It has become standard practice in lattice gauge theory to place the fields on a simple hyper- 
cubic lattice. While convenient for the purpose of doing numerical calculations, the restriction 
to hypercubic lattices is rather narrow and special from the perspective of continuum field 
theory. Since our interest ultimately is in defining and taking a continuum limit, we will set 
up our formalism on lattices more general than the hypercubic one. This prevents us of from 
making short cuts and helps us direct our attention to the proper structures. 

We take discrete space-time to be some ci-dimensional complex A built from oriented /c-cells. 
For the cases k = 0, 1, 2 these are referred to as sites, links (with a direction), and plaquettes 
(with a sense of circulation). It will often be convenient to view the /c-cells of A as the 
generators of abelian groups Cfc(A). Their elements, called /c-chains, are linear combinations 
of fc-cells with coefficients in Z. If c is a fc-chain, then so is — c; we get the latter from the 
former by reversing the orientation for all its &-ceils. For any reasonable choice of A there 
exists a boundary operator d, which is a linear operator d : Cfc(A) — > C/t_i(A) with the 
property dod = (the boundary of a boundary always vanishes). For example, the boundary 
of an oriented link 1 that begins on site nj and ends on site nj is dl = rif — m; which is 
the chain consisting of the 0-cells rif and nj with coefficient +1 and —1 respectively. The 
boundary of an oriented plaquette p is the chain of 1-cells lj surrounding it: dp = J2i 
where the plus/minus sign is chosen when the orientations agree/disagree. 

The following elaboration on language might help prevent confusion later on: each /c-cell of 
A comes with exactly one of two possible orientations; i.e. the choice of orientation for the 
fc-cells is fixed, albeit arbitrary. This means in particular that if the oriented link 1 is a 1-cell, 
then —1 is not a 1-cell (although it still is an oriented link). 

On a d-dimensional lattice or cell complex of this general kind, we are going to consider a 
gauge theory with partition function 



where the action functionals 5b and Sf will be specified shortly. As usual, matrices U taking 
values in a compact gauge group G are placed on the links of the lattice. More precisely 
speaking, a lattice gauge field configuration is a mapping from the 1— cells of A into G. The 
mapping extends to all oriented links by the convention U{—\) = f7(l) , which is motivated 
by the interpretation of the Vs as discrete approximations to the path-ordered line integrals 
of an underlying gauge field A: U(l) ~ Pexp(Jjyl). The a priori statistical weight of the 
lattice theory is a product of Haar measures dU over the 1-cells of A: [dU] = Y\idU(l). 
Concrete calculations will be carried out for the unitary groups G = \J(N C ). 
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In addition to the gauge fields, the theory has a second ingredient: complex bosonic and/or 
fermionic fields that are placed on the sites of A and transform according to the fundamental 
vector representation of the gauge group. The fermions are denoted by ip, the bosons by (p. 
They are vectors not only in iV c -dimensional color space, but also in a "flavor" space with 
dimension iVf (fermions) and N\, (bosons). 

These bosons and fermions are auxiliary (i.e. unphysical) degrees of freedom introduced solely 
for the purpose of inducing an effective action for the gauge fields. Unlike the conventional 
matter fields of lattice gauge theory, they do not propagate all over the lattice. Rather, each 
one of them is constrained to hop (in the presence of the lattice gauge field U) along the 
boundary chain of some plaquette. In other words, there is a one-to-one correspondence 
between the 2-cells p of A and sets of complex boson and fermion variables, {<p p } and {V'p}- 
The orientation of p determines the sense of circulation of the hopping of the auxiliary parti- 
cles {<£>p} and {^ p }; see Figure^ The integration measure [dip] [dtp] is taken to be the product 
of flat measures. 

To fix the precise details, we must distinguish between gauge groups of two types: those where 
the vector (U) and covector (C/ _1 ) representations are related by an inner automorphism, 
and others where they are not. The groups SO(N c ), 0(N C ) and Sp(2A r c ) belong to the former 
type, the unitary groups \J(N C ) and SU(A r c ) to the latter. For gauge groups of the former type, 
the setup we have described — one set of variables per 2-cell — would already suffice. For the 
latter type, however, and especially for \J(N C ), which the present paper focuses on, we have to 
double the set of auxiliary variables: it will be seen that, to induce a good effective action for 
the gauge fields, hopping must take place in both the clockwise and the counterclockwise sense 
for all plaquettes. Thus for each 2-cell p of A, we introduce two sets of auxiliary variables, 
one associated to the 2-cell with its proper orientation (+p), and another one where the 
orientation is reversed (— p). 

These words are put in formulas as follows. Let p = {p, — p}, and let L p denote the "length" 
of ±<9p, i.e. the number of links contained in the boundary chain of ±p. Then, fixing some 
oriented plaquette p, we write its boundary as a formal sum of L p oriented links with positive 
coefficients: dp = l ni ,n 2 +^n.2,n3 + ■■■■, where dl njtrij+1 = nj + i —rij. The 0-cells n\, . . . , ul p are 
the sites visited by dp, arranged in ascending order as prescribed by the sense of circulation 
of p. To take notational advantage of the cyclic structure of the boundary chain dp, we 
identify n\ = til v+ i. With these conventions, we put U p (rij + i, nj) = U(l njinj+1 ), and define 
the actions Sf and 5b to be 

S { [ip,ij,U] = ^2^2 ( m f.p^p( n i)^p(%) ~ $p( n i+i)^p( n i+i>%OV>p( n i)) > ( 2 ) 

±p j=l 

Lp 

Sb[p,<p,U] = ^2^2 ( m b,p£ P (%)<£ p (nj) - <p p (n j+1 )Up(nj +1 ,nj)<p p (njf) . (3) 
±p i=i 

To keep the expressions transparent, we have suppressed the color {i = 1,2, . . . , N c ) and 
flavor indices (a = 1, 2, ... , iVf /)-,); it should be clear how to restore them. For example, (pUip 

is short-hand for X^ C j 2 =i *l2a=i (p ll ' a U ni2 ip l2A . The notation X)± P m eans that each plaquette 
occurs twice in the sum, once each for the two possible orientations. The parameters mf iP and 
m^ p are referred to as the (local) fermion and boson masses. We allow for the possibility that 
they depend on the plaquette label p in general, but on a lattice with translational invariance 
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Figure 1: Hopping of the auxiliary bosons and their coupling to the gauge field. 

and no external gravitational field, we will take them to be constant. The field integral makes 
sense for any set of m^ p G C, while the local boson masses must satisfy Remt, iP > 1 for 
convergence. 

The lattice theory so defined is intended to be a discretization of (i-dimensional (Euclidean) 
Yang-Mills theory with gauge group G = U(iV c ) or SU(-/V C ). In the second half of the paper 
we will eliminate the gauge field so as to pass to a dual description. For the time being we 
address the question whether (JJ) has a critical point with a continuum limit equivalent to 
Yang-Mills theory. 

The field integrals for the auxiliary fermions ip, tp and bosons (p, <p are Gaussian. Carrying 
them out we obtain a product of determinants (for fermions) and inverse determinants (for 
bosons), one each for every oriented plaquette p. The determinants from the fermions asso- 
ciated with p are Det^ (m f £ — U(dp)), with U(dp) being the ordered product of the U's 



The corresponding factor from the same plaquette with the opposite orientation, — p, is just 
the complex conjugate of this. The Gaussian boson integrals give a similar answer except 
that the determinants in that case go in the denominator. Thus, combining factors we have 



Note that this is a product over 2-cells, i.e. each oriented plaquette p occurs only once (and 
the result does not depend on the orientations chosen). 

By sending the result of the integration back to the exponent and dropping an irrelevant 
constant, we obtain 



along the boundary chain of p: 



U(dp) = U p (ni,nL p ) U p (n Lp , n Lp _i) • • • £7 p (n 3 ,n 2 ) U p {n 2 ,ni) . 





(4) 
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with fermion and boson induced actions 

S- md , f [U] = -27V f Re^Trln(l-Q fiP C/(a P )) , (5) 
p 

Sind,b[U] = +2N h Re ]T Tr In (l - a b , p U(8p)) . (6) 
p 

In place of the fermion and boson masses, we have introduced the coupling parameters «f p = 

m { p p and a^ tP = m b p . We take the range of these parameters to be — oo < af tP < oo and 
— 1 < a^, p < 1. Note that the partition function © becomes singular in the limit a^ iP — > 1, 
since the function U h-> |Det(l - aU)\~ 2Nh does so on the codimension one submanifold of 
matrices U with at least one eigenvalue equal to unity. 

There is an obvious way |0j in which to approach Wilson's lattice gauge theory |2.S] 1241 125] 
with action 

s w [^] = -|-E ReTr ^p) (7) 

iVc p 

from © and ©: specializing to the case of a hypercubic lattice, we take the couplings a^ ;P 
and af iP to be independent of p. We then send them to zero, and the number of flavors iVf/b 
to infinity, while keeping the products af/b-Wf/b fixed. Since limAr^oo iVln(l + x/N) = x, this 
limit exactly reproduces the Wilson action (jjj), with /3/iV c = 2(A r b a b — N{(X{). If only fermions 
are used, the coupling «f = m f p must be negative, corresponding to a complex mass, say 
mf = |mf | e 17r / Lp . 

It is very unclear to us, however, whether this double limit of an infinite number of infinitely 
heavy particles is going to lead to any advance in our knowledge about Yang-Mills theory. 
Certainly, the case of a single boson flavor (iV b = 1) treated in 0| bears no relation with 
Yang-Mills theory, let alone QCD. 

In the present paper we are going to propose and investigate a more interesting possibility. 
The partition function (JIJ becomes singular when all a^ jP are sent to unity (say, uniformly 
in p), which is in fact a critical point with a diverging correlation length. Our main message 
will be that, on suitable lattices and in high enough space-time dimension, the critical theory 
with a bjP — ► 1 (and Of jP ^ 1) is expected to be in the universality class of U(-/V c ) Yang- 
Mills theory, provided that iV b > N c . The universality conjecture can readily be checked in 
dimension d = 1 + 1. There, it will be demonstrated that universality holds for iV b > N c , 
while it fails for Ny, = N c . 

3 The Critical Point = 1 

3.1 Lattice versus continuum 

To communicate our argument, we must first recall the general scenario [241 125] by which 
lattice gauge and continuum gauge theories are related. Imagine some smooth ci-dimensional 
configuration of the gauge field A, and superimpose on it a fine mesh in the form of some 
d-dimensional lattice or cell complex A with oriented elementary plaquettes p. The rule of 
translation from the continuum gauge field A to the lattice field U is given by 

U(dp) = Pexp <b A = l + F(p) + ... , (8) 

Jdp 



6 



where P means path ordering, and F(p) is the field strength (the curvature of the gauge 
connection A) evaluated on p (in the standard setting of a hypercubic lattice, this means that 
if p is a plaquette based at the site n and is parallel to the /i^-plane, then -F(p) = a^F^^n) 
with a the lattice constant). 

Clearly, the lattice approximation will be reasonable if the mesh formed by the lattice is 
fine enough in order for field curvature effects to be small on the lattice scale. Under such 
conditions P exp ^ A, and hence U(dp), will be close to unity for every elementary plaquette 
p. Thus, if the lattice approximates the continuum, the plaquette matrices U(dp) fluctuate 
only weakly around the zero-field strength configuration U(dp) = 1. 

Conversely, if the statistical weight of the lattice theory sharply peaks at unity U (dp) = 1 for 
all p, then the lattice theory has a very large correlation length and is close to a continuum 
limit. We may then pass to a continuum field theory formulated in terms of the field strength 
F, by expanding around U(dp) = 1 and using the correspondence (JHJ). 

The main principle then is this. Let the lattice gauge theory be given by a product statistical 
measure Y[ p w t\U(dp))[dU], where t is some coupling parameter with critical value t c . We 
require the weight function wt(U) for t close to t c to be very strongly localized at the unit 
element so as to make excursions away from unity statistically very rare. More precisely, 
we want wt c (U) = S(U) to be the Dirac <5-function supported at the unit element of the 
gauge group, and wt(U) for t ^ t c to be some smeared version thereof. When t moves to t c , 
the smearing is undone, the lattice statistical weight approaches the product of ^-functions 
Y[ p 5(U(dp)) , the correlation length goes to infinity, and the lattice gauge theory converges 
to a continuum limit. Under favorable conditions, this continuum limit will be Yang-Mills 
theory with action functional — Sym oc J Tr F A *F. 
This scenario is realized for the Wilson weight function 

when the parameter f3 = t is sent to t c = oo. For one thing, one easily shows 

lim I f(U)wJ(U)dU = f(l) 

U(iV c ) 

for any smooth function / on U(iV c ) (or another compact gauge group, for that matter), 
which is equivalent to saying lim^oo vj^ (U) = 5(U). Moreover, it is strongly suggested by 
numerical simulations that the continuum limit approached by Wilson's lattice gauge theory 
for (3 — > oo has all the properties expected of quantum Yang-Mills theory in the continuum. 
The ultraviolet stability of Wilson's theory has been established by rigorous analysis |26| . 

There is of course nothing unique about the Wilson action and many other one-parameter 
families of weight functions do the same or an even better job. Let us single out one example. 
In the representation theory of compact semisimple Lie groups G there exists a statement, 
called the Peter- Weyl theorem |27j . which implies that the Dirac <5-function on G can be 
built up from the complete set of irreducible representations D A of G as follows: 

5{U) = vol(G)" 1 Y, d x TrD x (U) , (9) 

all A 
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where d\ = TrD A (l) is the dimension of the representation. Hence, if 02(A) > is the 
quadratic Casimir invariant evaluated in the representation D x , the function 

wf K {U) = vol(G)" 1 ^ e~ C2[X)t d\ TvD x (U) (t > 0) (10) 

all A 

approaches the 5-function as t — > t c = 0. For this reason one expects it to provide (in the 
limit of small t) a valid lattice regularization of Yang-Mills theory with gauge group G. 

This weight function is called the heat kernel: it solves the heat equation dtWt(U) = Awt(U) 
(A being the Laplace-Beltrami operator) on G with initial condition lim^o+ u>t{U) = <5(£7). 
The "time" parameter t is the analog of 1/(3 for the Wilson weight function. In the abelian 
case G = U(l), the model with weight function is known in condensed matter theory as 
the Villain model. 



3.2 5— function limit 

After these preparations, we are ready to get to the point: we are going to investigate the 
determinantal weight function 

w Nbta (U) = \Det(l-nU)\~ 2Nb / J |Det(l - a U)\~ 2Nh dU (11) 

of the boson induced gauge model © with ctf = and a = < 1 close to unity. We will 
be interested mostly in the case G = U(iV c ), but note that the definition of wjy h:a still makes 
sense if G = U(iV c ) is replaced by any of its closed subgroups. 

As a warm up, we will look at two special cases. The simplest example is G = U(l) with 
iVt, = 1, U = e 10 and Haar measure dU = dO. In that case, elementary manipulations show 

l-ae 1 " 2tt ^-^ 

1 1 nez 

from which it is seen that w 1 Q approaches the Dirac 5-function on U(l) in the limit a — > 1. 
The difference to the Villain model (jlUj) is that the Gaussian cutoff e _C2 ^ n ^ = e - ™ 2 * has been 
replaced by an exponential cutoff e - '™'*, with t = ln(l/a) > 0. 

Another nice example is G = SU(2), still with N\, = 1. Again, by straightforward manipu- 
lations on the determinantal weie ht function JUJ) evaluated on U = e i0CT3 G SU(2), one finds 

with /su(2) dU = 27f2; 

Wl ' a[G >-V*> |l- a e i9 | 4 " 2^^ a ln + ij sine 

1 1 n=0 

On the right-hand side we recognize the character xs( el(T3 ) — si n ((25" + 1)8) / sin 8 of the spin 
5 = n/2 representation of SU(2), and its dimension 2S+1 = n+1. Both sides extend uniquely 
to class functions of SU(2), i.e. to functions on SU(2) that are invariant under conjugation 
U 1 ^ gUg" 1 by g £ SU(2). As a result, we deduce lim a _»i wi >a (U) = 5 SV ^(U) from ©. 

Fact. Let f : U(iV c ) —> C be an analytic function. Then, for any N\, > N c , 

lim / f(U)w Nb , a (U)dU = f(l) . (12) 
U(JV C ) 
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Thus we are claiming the desired ^-function property for > N c . The proof will be given 
in the next subsection. For now we make two comments: (i) There is nothing special about 
\J(N C ) in this context, and we expect a similar statement to be true for each of the compact 
matrix groups G = SU(iV c ), SO(N c ), 0(N C ), Sp(2iV c ). (ii) In view of the Peter-Weyl theorem, 
the statement Q12JI implies that the complete set of irreducible representations of U(iV c ) occur 
in the character expansion of WN bta (U) for N\, > N c . Conversely, one can show that some 
representations are missing for N\, < N c , which means that the inequality N\, > N c in the 
statement cannot be relaxed but is optimal. 



3.3 Proof of fact 

We will give an elementary proof, and start with a few considerations that simplify it. 

First of all, if the function U i— > |Det(l — a U)\~ 2N is concentrated near unity U = 1 for some 
value of the exponent N, then it will be even more so for exponents larger than N. Thus, 
if the statement is true for some value N of N^, it will certainly be true for all values of 
greater than N. It therefore suffices to show that the statement holds for the border line case 
iVb = N c . We abbreviate the notation by writing w a (U) = WN c ,a(U). 

Second, since the weight function w a is invariant under conjugation U i— > gllg -1 , the operation 
of replacing the analytic function / by its average / av over conjugacy classes, 

r(U)=vol(XJ(N c ))- 1 J f(gUg- l )dg, 

U(7Vc) 



does not change the value of the integral (|12[). We may therefore assume / to be invariant. 

Third, given invariance under conjugation, we may view w a and / as functions on the maximal 
torus T = U(l) JVc parameterized by the eigenvalues e , . . . ,e l6,JVc of U, and we may reduce 
the integral over U(iV c ) to an integral over T. Let J(e w \. . . ,e i6,JV <=) = l\ k<l \e i9k - e Wl \ 2 be 
the Jacobian of the polar coordinate map (U(-/V c )/T) xT-t \J(N C ). Then we have 

(f)a 



f(U)w a (U)dU 

U(N C ) 



<1>« 



where the angular brackets mean 



(f)a= I f(e ie \...,e l6 ^) ^ ) —d0 1 ---dB Ne 

Tl -i I 1 ~ ae A 2Nc 



l<Vc 



[0,2vr; 

To prove the fact (|12(1 we must show 

lim-^ = /(l,...,l) (13) 

for all analytic functions / on the maximal torus T = \J(l) Nc which extend to functions on 
U(iV c ). (Such functions on T are invariant under the Weyl group of U(iV c ), i.e. they do not 
change under permutations of their arguments.) 

We will actually establish the limit (|1 3|) for the larger class of all analytic functions F : T — > C. 
By definition, such functions are absolutely convergent sums of the basic functions e 1 ^ nfc6 ' fe 
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with integer exponents n&. It therefore suffices to prove 1)13(1 for the complete set of these 
basic functions. So let 

J? = e i(ni9i+...+njv c 6jv c ) 

with any (n\ , . . . , njv c ) £ Z c . Without loss we may assume that the ordering of variables has 
been adjusted so that the first p integers m, . . . , n p are positive or zero, while the last N c — p 
integers n p+ i, . . . , n/v c are negative. 

Now we evaluate (F) a in the limit a — > 1. The first step is to switch to the variables z k = e l9k , 
which yields the expression 



m n Nc 

.Gl * * * /C 



(F) a = (i/a) N - <b - ^ T\(z k - zi) 2 d Zl 

J Y\f=i( z j-^ Nc ( z j-»- l ) Nc 



U(l) 



A'c 



■ ■ ■ dz 



By the signs assumed for the integers m , . . . , njv c , the integrand is obviously regular at zero in 
the complex plane for each of the variables z\, . . . ,z p . Simple power counting shows that the 
same is true at infinity for z p+ i, . . . , zn c - (This is to say that the integrand for k = p+1, . . . ,N C 
decays as z^ 2 or faster at infinity.) Therefore, our strategy now is to contract the contour 
of integration to zero for the first p variables, and expand it to infinity for the last N c — p 
variables. In doing so, we pick up contributions from the poles of order N c at z\ = . . . = z p = a 
inside the unit circle U(l) C C, and at z p+ \ = . . . = zn c = a -1 outside the unit circle. By 
the residue theorem, we then arrive at the exact formula 

(F) a = (i/a)^(2^(-27ri)^-f(iV c -l)!-^ 

*N C llk<l 



2 p+i = 



This expression is divergent at a = 1. The highest-order singularity, (a — a -1 )"^, occurs 
when all derivatives act on nfc</( z fc — z i) 2 or on the denominator. The order of the singularity 
is reduced if one or several of the derivatives act on the monomial F = z™ 1 • • • z N Nc . To 
compute the limit a — > 1, it is enough to retain the leading-order singularity. The leading- 
order singularity is picked up by evaluating F at z\ = . . . = z p = a and z p+ \ = . . . = z^ c = 
a" 1 and taking it outside of the expression. What is left behind is just the value of the 
integral {l) a obtained by replacing F by unity. Hence 

(F) a = aJ ni|+ "- +|njv cl (l) Q + less singular terms . 
So we conclude lim Q _ + i(F) Q /(l) a = F\ _ _ _ , and the proof of ()12() is complete. 



4 Continuum Limit in Two Dimensions 

The result (fT2|) ensures that on sending all coupling parameters ab, p — > 1, the induced gauge 
model © for > A^ c becomes critical, which allows a continuum limit to be taken on any 
reasonable direct system of lattices. Based on universality, we expect this continuum limit to 
be quantum Yang-Mills theory, at least generically. 

A precise investigation of the universality conjecture can be made, and will now be made 
using harmonic analysis on the gauge group, in the simple case of d = 1 + 1 dimensions. 
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There, the universality mechanism at work is the central limit principle in its basic form: 
computing the two-dimensional theory essentially amounts to taking convolutions of the 
plaquette distribution WN h) a{U)dU with itself and, by a central limit theorem for compact 
Lie groups, multiple convolution sends a large class of distributions to the universal heat 
kernel family w^ K (U)dU written down in 

In short, the central limit principle we will exploit is this. (Although we are going to pursue 
the case of U(A C ), we here state the principle for a semisimple compact Lie group G. The 
extension to U(A C ) will cause minor complications.) Let wt(U)dU (t > 0) be a one-parameter 
family of smooth AdG-invariant distributions on G such that limt_>o w t(U) = 5(U). Using the 
exponential mapping Xi->[/ = exp X we can pull back the family to a family of distributions 
dfit(X) on the Lie algebra of G (or, rather, to a domain of injectivity of exp in LieG). With 
respect to d(it(X) we compute the expectation of the Killing form (X, X) = — Tr ad(X)ad(X). 
If there exists a "diffusive scaling" limit, i.e. the expectation of (X,X)/t stays finite when t 
is sent to zero, then a central limit principle applies: denoting the iV th convolution of wt with 
itself by w* N , we have Huin^ooW*^ = w^ K . 

We will see that the plaquette distribution of the boson induced model © for Ny, > N c + 1 
satisfies the diffusive scaling criterion, with t ~ (1 — «b,p) 2 - This will allow us to take a 
continuum limit which can be considered as a rigorous definition of 2d quantum Yang-Mills 
theory by the reasoning of Witten |28j . 

On the other hand, for = N c the diffusive scaling criterion turns out to be violated! A 
continuum limit can still be defined, owing to (|12|). This, however, is not Yang-Mills theory 
but an unusual theory which, in the "first-order formalism" with an auxiliary LieU(A c )- 
valued scalar field 4>, is given by an action functional 

S = -i f Ti(f)F + n [ || <j) Id d 2 x , (14) 
Jt, Jt, 

where || <j) \\\= Y2f=i I'Aj'Ij the 4>j being the eigenvalues of <j). This exists as a renormalizable 
theory because the Cauchy distribution on u(N c ) = LieU(A" c ) approaches under subdivision 
a distribution which is stable, and yet different from the heat kernel. 

We shall begin to substantiate these assertions in Section 14.21 



4.1 One-Plaquette Model 

Before we undertake the study of the two-dimensional models where Yang-Mills universality 
rules, we dispose of those where it does not. As a simple test, we look at a cell complex 
consisting of a single plaquette p and consider the expectation value of Tr U = TvU (dp): 

Z(a f ,a h ) J |Det(l - a h U)\ 2N b 

U(7Vc) 

In U(iV c ) Yang-Mills theory the expectation of TrC/(C) (the holonomy along any loop C) 
goes to Tr 1 = N c when the coupling is sent to zero. The same happens with VF(af,ab) 
in the limit «b — > 1, Ab > N c , for in that case the statement (|12|) applies, and gives with 
f(U) = TrC/|Det(l - a{U)\ 2Nt : 



lim W(a { ,a h ) = TrC/| = A c (N h > N c , a { / 1) . 
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While this is a statement made for the fundamental representation, (|12j) asserts that a similar 
result holds true for the trace Tr D X (U) in any representation D A . 

What happens in the other cases? We separately look at the boson induced models (af = 0) 
with N\> < N c , and at the fermion induced models («b = 0), starting with the former. Using 
complex contour integration and residue calculus, we show in App endix I A . 21 that the following 
holds true: 

lim W(0,a h ) = N h (N h < N c ) . 

Thus in that case W(0, 1) falls short of reaching the maximal value N c allowed by the bound 
|Tr£7| < A^. The general dependence can be computed numerically (see Appendix IA.2|) . and 
is shown in Figure [21 for the case N c = 3. An interesting observation, proved in Appendix 
fA~2l is that for N h < N c the function W(0, ab) is exactly linear: W(0, «b) = A^a^. 

Turning to the fermion induced case, from the invariance of dU under £/—»—{/ and U — > U , 
we deduce W(a{,0) = — W(— «f,0) = W(l/af,Q), so it suffices to restrict attention to the 
interval — 1 < at < 0. Figure |3 shows numerical results (Appendix IA.1|) in that range. We 
observe that W(cef,0) is a monotonically increasing function of — Of for |cuf | < 1. The global 
maximum, attained at ai = —1, can be computed analytically by a fermionic variant of the 
method of Howe pairs described in Section 14.51 

m-i.o) *»• 



N, + N c 

Thus, the maximum again falls short of reaching the limit posed by the bound |TrZ7| < A^ c . 

These results show that in the fermion induced model with any finite number Nf of species, 
and in the boson induced model with A^ < A^ , the single-plaquette action cannot ever enforce 
complete suppression of the fluctuations of U away from the unit element. Consequently, we 
do not expect a continuum limit of Yang-Mills type in these models. 



4.2 Character expansion 

We now embark on a detailed study of the boson induced gauge model © with at least as 
many flavors as colors (A" c > A^). As a preparatory step, we do some harmonic analysis for 
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Figure 3: Same as Figure |21 but now for the fermion induced model. 

the one-plaquette distribution 

dfi Nhia (U) = |Det (1 - aU) \~ 2Nb dU . (15) 

Our motivation is that we wish to compute multiple convolution integrals of this distribution 
with itself, and transforming to the appropriate Fourier (or harmonic) space turns convolu- 
tions into simple multiplications. 

Recall that by the Peter-Weyl theorem any L 2 ~function on a compact Lie group can be 
expanded in the matrix entries of its complete set of irreducible representations D x . Since 
the distribution function in (|15[) is a class function of U(iV c ), the expansion proceeds by U(iV c ) 
characters X\(U) = Tr D X (U): 



|Det(l-aC7)| 2Nh = c x (a) X x(U) , (16) 



-2N b = 

all A 

and, by the orthonormality of characters, the expansion coefficients are 



c x (a)= J \Det{l-aU)\ h x\{U~ l )dU . (17) 

U(iV c ) 

Here, and throughout this section, we normalize Haar measures by Jjj, N \ dU = 1. From 
statement (|12j) we already know the limit 

lim i — TT = Xa(1) = <^A , 
l c (a) 



a 



which asserts convergence of the ratio of coefficients to the dimension of the representation A. 
The goal of the current subsection is to gain a more precise understanding of the expansion 
coefficients c\(a) close to a = 1. 

To that end we will use the Cayley map 

' y ' 1-X 
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to express the integral (|17|) over U(iV c ) as an integral over the anti-Hermitian N c xN c matrices, 
X £ u(iV c ). The Jacobian of the Cayley map is ~Det~ Nc (l — X 2 ), which is to say that there 
exists a flat positive density dX on u(N c ) such that 

7*(<fl7) = Det _JVc (l - X 2 ) dX . 

Note —X 2 > for X £ u(N c ). We then immediately have the statements 

J Det~ Nc (l -X 2 )dX = j (#7=1, 

u(7Vc) U(7V C ) 

J Det _JVb (l - X 2 )dX < 1 for N h > N c , and (18) 

u(7Vc) 

J Tr(-X 2 ) Det _Arb (l - X 2 ) dX < 1 for N h > N c + 1 . (19) 

u(7Vc) 

The last of these follows from the elementary inequality Tr(— X 2 ) < Det(l — X 2 ). 
Making two variable transformations in sequence, 

u = \±±, x = ^Xeu(JV c ), 

1 — x 1 + a 

we bring the integral formula for the expansion coefficients into the form 

N? 



c x (a) = (l- a) - 2 ^( 1 -^) C I Xx((l-x)/d + »)) 
M 7 V ^ \l + aj J Det^-^ (1 - x 2 ) 

u(N c ) 



Det- Nh (l -X 2 )dX 



In view of the finiteness statement (|18)) it is clear that the integral for c\(a) will localize at 
unity U = 1 in the limit a — > 1, for JVb > A^c. Hence, to get an accurate approximation to 
these c\(a) near a = 1, the only further input we need is an understanding of the characters 
X\(U) close to unity. 

Let A denote a set of N c integers, A = (Ai, A 2 , . . . , Ajv g ), let e ie = (e Wl ,. . . , e i0Ar -) G \J(l) Nc , 
and define the elementary antisymmetric torus function £\ by 

£ A (e ie ) = sgn(7r)e iAie -( 1 ) e iA2e -( 2 ) ••• e Uw « e *( w «) , 

i"65'jv c 

where the sum is over all permutations of {1,2,... ,A' r c }. By a classic result of Weyl |29j . 
the irreducible representations of XJ(N C ) are in one-to-one correspondence with ordered sets 
A, Ai > A2 > • • • > Atv c , and the character associated with A is the class function xx(U) 
determined by 

Xx(e w ) = Cx +P (e ie )/C P (e W ) (20) 

where p = {N c — 1, N c — 2, . . . , 1, 0). The character at unity computes the dimension of the 
representation, which by the Weyl dimension formula is 

Xx(e°) = d x = A(X + p)/A(p), (21) 
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^■W = nfc<i(-^fc — A/) being the Vandermonde determinant. 

Because characters are joint eigenfunctions of the full ring of invariant differential operators, 
the characters x\{U) of U(iV c ) = U(l) x SU(iV c ) with central factor U(l) separate. If we put 
U = e x , a U(l) factor e(i/ N <:) TTX with charge q(X) = Xj splits off. Moreover, since the 
SU(iV c ) part only depends on the traceless part of X £ u(N c ), the characters expand around 
unity as 

( X\ 

I = e (q/N c )TrX x fSV(N c ) , Q X . e -TrX/N c j 

with some coefficient A = A(X). In parallel, one can directly expand Weyl's formula (|2()|) 
for x\ t° second order in the angles 6j. By doing so and comparing coefficients (or by more 
sophisticated techniques not discussed here), one finds 

A(X) = (N 2 - I)" 1 (Cas 2 (A) - q(X) 2 /N c ) , 

where Cas2(A) is a quadratic Casimir element of U(iV c ), 

Cas 2 (A) = x j( x j +N C + 1- 2j) , (22) 
j"=i 

which is associated with the invariant quadratic form — Tr(X 2 ) in the standard way |27j . 
Having gathered enough information about the characters x\ near unity, we now return to 
the task of computing the behavior of the expansion coefficients c\(a) close to a = 1. If we 
insert the small-X expansion of x\( eX ) the integral formula for c\(a) and take the ratio 
c\(a)/co(a), we encounter the second moments of the distribution Det~ Arb (l — X 2 )dX: 



E Nh Tr(X' 2 ) = Z- 1 J Tr(A' 2 )Det- JVb (l - X 2 ) dX , (23) 

u(JV c ) 

E Nb (TiX) 2 = Z- 1 J (TrX) 2 Det _Arb (l - X 2 ) dX , (24) 

u(iV c ) 

with Z = f, N , Det _JVb (l - X 2 ) dX. These are finite for N h > N c + 1, by the inequality 
(UU) and < -(TtX) 2 < -N c Tr(X 2 ). The first moment E Nh TrX vanishes by parity. It is 
therefore easy to prove the following statement: 



d x (l-^(l-a) 2 (B iq (X) 2 + B 2 C™ 2 (X))+K(a)\ , (25) 



Fact. The coefficients c\(a) for > N c + 1 have a Taylor expansion 

c\(a) _ j ( 1, 
c (a) 

with a remainder term TZ(a) that vanishes faster than (1 — a) 2 in the limit a — ► 1. The 
coefficients B\ and B 2 are determined by the linear system of equations 

-E Nb (TvX) 2 /N c ^ = (N C l\ (26) 



E,v b Tv(X 2 )/N c J V 1 N C J \B 2 



The leading singularity in the normalization cq{q) is 

c (a) ~ (1 - a)~ 2N ^ +N * . (27) 
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4.3 Continuum limit 

We now place the boson induced gauge model © with N\> > JV C + 1 on a two-dimensional 
cell complex A approximating an orientable compact Riemann surface S. Such a complex 
consists of plaquettes, links and sites, with every link joining at most two plaquettes, and the 
plaquettes receiving an orientation from S. 

To each plaquette p of A, we assign an area A p determined by some choice of Riemannian 
structure of S. Measuring area in units of a fundamental area a 2 , we put a p = 1 — ^jA~ p ja. 
This specifies the set {a} of coupling parameters of our boson induced gauge model on A. 

We now focus on the partition function 

Z A ({«})= [[dU]H\Det(l-a p U(dp))\- 2Nb . 
J p 

The goal is to demonstrate that, when the cell complex A is refined so as to approximate X 
ever more closely, Z\({a}) converges to the known partition function of XJ(N C ) Yang-Mills 
theory on E, with a particular choice for the U(l) coupling. 

Since any £ can be made by gluing of surfaces with disk topology, the first step we take is 
to compute the (boundary-value) partition function, T\(U,{a}), for a two-dimensional cell 
complex A approximating a disk, where 

U = U(dA) = ]J U(l) 
ledA 

is the product of U's along the boundary <9A, and these boundary £/'s are held fixed to 
prescribed values. To compute T\(U,{a}), we need to integrate over all the matrices U 
associated with the links in the interior of A. These link integrations can be carried out 
in any order, and because A has disk topology (i.e. every interior link joins exactly two 
plaquettes), it suffices to show how to do one of them. 

Hence, let 1 be any interior link of A, and let p and p' be the two oriented plaquettes joined 
by 1. We are going to do the integral over the link matrix U(l). For that purpose, let the 
holonomy along the boundary chain of p be 

U (dp) = U p (m,n Lp ) ■ ■ ■ U p (n 3 , n 2 ) U p (n 2 , ni) , 

and similar for p', where our notational conventions are as specified in Section |2j Without 
loss we may assume that 1 is the link joining the sites ri\ = n' 2 and n 2 = n'<y. 

U = U(\) = U p (n 2 , ni ) = {U^in^n'^y 1 . 

Abbreviating the notation by putting 

[/(dp) = WU , U{dp') = VU~ l , 

we are then faced with the convolution integral 

I(V,W)= J |Det(l - cyW/" 1 )!" 27 ^ |Det(l - u p UW)\~ 2Nh dU . 

U(7Vc) 
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To carry it out, we use the character expansion |Det(l — aU)\ 2Nb = ^2\C\(a)x\{U) and 
the fact that characters reproduce under convolution, 

Xx(VU~ 1 ) X y(UW)dU = 5x, x , XxiyW) , (28) 

U(JVc) 

to obtain 

A aX 

We now iterate the procedure, and successively do all inner link integrations using the convo- 
lution law (|28|) for the characters. The resulting expression for the boundary- value partition 
function is 

where U = U(dA) = J7a(^l A) til a -i) • • • U\(n2,ni) is the holonomy along <9A, as before. 

Here is the point where we take the continuum limit. If we refine the lattice discretization, 
and keep on refining it so that A becomes a closer and ever closer approximation to a disk 
D (or some other surface D diffeomorphic to a disk), the elementary plaquette areas A p go 
to zero, the parameters a p = 1 — y^A^/a approach unity, and we may eventually use the 
asymptotic law (|25|). giving 

/ A \ 

co(a p ) + . . . , 



n ^ = n (i - (*i ^(^) 2 + b * 



with corrections that vanish in the limit. Since limAr_ >00 (l — x/N) N = e~ x , the plaquette 
areas A p exponentiate and piece together to yield the total area ^ A p of the surface D. To 
eliminate one irrelevant constant from our expressions, we adopt the convention of measuring 
area in suitable units: fi = (B2/CL 2 ) x ^2 p A p . Taking the continuum limit then leads to 

T d (U,im) = ^2^ ex P ("f (Cas 2 (A) + (B l /B 2 )q(\f)) Xx {U) . (29) 

A 

We have dropped a (diverging) multiplicative constant np c o( a p) that arose from our using 
an unnormalized statistical weight function. (Physically speaking we made use of our freedom 
to set the vacuum energy of flat space to zero.) 

An expression for Td(U, fx) of the given form constituted the starting point of Witten's com- 
binatorial treatment From it, he computed (among many other things) the Yang-Mills 
partition function Z g (fi) for any orientable compact Riemann surface of genus g and dimen- 
sionless area /i. Instead of repeating that computation here, we just quote the answer it 
entails in the present case: 

Z g (p) = Y / d- 29+2 e-^W+( B ^ B ^ 2 )/ 2 , (30) 

A 

and verify it for the two simplest examples, the sphere and the torus. 

A sphere (g = 0) with area \i can be made by gluing together two disks D + and D_, say 
with areas //+ and fx- = fi — The boundary-value partition functions Td + and 
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depend only on the total holonomies along dD + and <9D_ respectively. If the two disks are 
to fit together to give an oriented sphere, we must have dD + = — <9D_, which means that the 
holonomies are inverse to each other: U{dD + ) = (/(dD )^ 1 . Thus, introducing U = U{dD + ) 
as the variable of final integration, we obtain 

Sphere (m) = J r D+ (u,[j, + )r D _(u~ 1 ,n-)du . 

U(JV C ) 

Inserting (|29j) and doing the resulting integral over a product of two characters with the help 
of ([28)1 . we indeed get the answer IjMOJI for genus g = 0. 

To manufacture a torus (g = 1) with area /x, we take a rectangle Q with the same area, and 
identify opposite edges. Writing the holonomy along dQ as a product over its four edges, 
U(dQ) = UVU" l V" 1 , we have 

Ztorus(/i)= J f Y Q (UVU- 1 V- 1 ,ii)dUdV . 

U(7Vc) U(7Vc) 

To do the first of these two integrals, we use 



/ xxiUAU^B) dU = X \(A)xx(B)/d x , 
Jg 



which is a consequence of the orthogonality relations obeyed by the matrix entries of an 
irreducible representation D x . Doing the second integral with (|28|). we again reproduce the 
formula (|30j) . now with genus g = 1. 

Let us summarize. From the (5-function property (|12[) we always expect the boson induced 
gauge model (JHJ) with iVb > N c to admit a continuum limit. In two dimensions — and on a 
direct system of two-dimensional cell complexes converging to a compact Riemann surface 
— , making the natural assignment (1 — a p ) 2 ~ A p and taking the continuum limit A p — > 
sends the disk boundary- value partition function of the boson induced gauge model to 
the expression (|29jh if iVb > N c + 1. By the reasoning of |28j . this expression is the proper 
local data to use for U(iV c ) Yang-Mills theory, and gives the partition function (|3U|) . 

More precisely speaking, the situation is this. The gauge group U(iV c ) is not semisimple; the 
number of its quadratic Casimir invariants is two (as opposed to one for a semisimple Lie 
group), and therefore XJ(N C ) Yang-Mills theory is not unique but comes as a one-parameter 
family. The procedure of canonical quantization associates the Casimir invariants Cas2(A) 
and q(X) 2 with the action densities Y.\ F %? d2x = — Tr(F A *F) and - £ Fj^FjUd 2 'x = 
— (TrF) A*(TrF), respectively. Thus the combination 

Cas 2 (A) + ( J B 1 / J B 2 )g(A) 2 

in (j30j) is the Hamiltonian arising from the action functional 

-S = Tr(F A *F) + ^- ^(TrF) A *(TrF) 

by canonical quantization. The Lagrangian depends on the parameters B\ and E>2, which in 
turn are given by the second moments of the distribution Det~ Nb (l — X 2 ) dX; see (|2fi|). 
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Note that this distribution becomes Gaussian for Ny, — ► oo. In that limit, — E7v b (TvX) 2 ~ 
N c /2N b and -E Nh Tr(X 2 ) ~ N 2 /2N h , and we get J3i/J3 2 = 0. Thus, the coupling Bi/S 2 
gives some measure of how much the distribution Det _Arb (l — X 2 ) dX differs from a Gaussian 
distribution. For iV c = 2 an easy computation gives 

-E Nh Tr(X 2 ) = — i— + — J— , -E Wb (TrX) 2 - 



2N h - 5 2iV b - 3 ' b v ; 2N h - 5 2iV b - 3 ' 

which yields Bi/JB 2 = |(iV b - 2) _1 . Thus in this case the ratio B 1 /B 2 for N h > N c + 1 = 3 
is always positive. 

4.4 The Cauchy case iVb = N c 

The continuum limit taken in the previous subsection does not go through for iV b = N c . The 
main obstacle is that the second moments of the Cauchy distribution Det _Arc (l — X 2 ) dX do 
not exist. The nonexistence can be verified by power counting, and is compatible with the 
failure of the inequality (|19|) in that case. It is also signalled by the expectation values of 
the positive quantities — Tr(X 2 ) and — (TrX) 2 we have just given for iV b > N c = 2, which 
are seen to become negative at iV b = N c = 2. Thus, although we still have a good integral 
formula, 

c x (a) = (1 - a 2 )~ N ? J X A fellf) Det_JVc C 1 " **) dX < 

we can no longer extract the behavior of c\(a) near a = 1 by expanding the argument of 
the character \\ under the integral sign, as this would lead to a divergent integral. Hence a 
different approach is needed. 

To get some inspiration, we turn to the simple case N c = iV b = 1 with distribution 



d9 



dx 



fo^ ll-ae*!-* - = (i7r) ~ 1(1 " ^^-t ( x E 1 



l + a + (l-a)x \ J \ ^ i Y — X 2 

1 + (1 — a)x 



The primitive characters of U(l) are Xn(e l9 ) = e mS (n £ Z). The coefficients c n (a) can easily 
be written down in closed form: 

c n ( a ) = (1-a 2 )" 1 ^™ 1 . 

Obviously, the effect of having a divergent second moment here is to cause a nonanalyticity 
in the Fourier variable n at n = 0. In the previous case of iVb > N c + 1 = 2, the ratio 
c n (a)/co(a) had zero derivative at a = 1 by parity symmetry (x \-> —x), but now 

d c n (a) 

lir ? w TT = n ' 

a-»l- aa co(aJ 

which forces a different scaling of the parameters a p , and thus a different continuum limit. 
We are going to show that the same situation occurs for all iV b = N c > 1. 

Consider the formula 

r A (a) = 44 = / XA f ^ft^ l Det"^ (l - X 2 ) dX , 
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which expresses the ratio r\(a) as an expectation of the character xx with respect to the 
normalized Cauchy distribution Det _Arc (l — X 2 ) dX. Because |xa(^0I is bounded from above 
by the dimension d\, the integral makes sense and defines r\(a) as a continuous function of 
a in the range — 1 < a < oo, and we have 

l r A(«)| < d\ ( — 1 < a < oo) . 

r\{a) is maximal at a = 1, where the upper bound r\(l) = Xa(1) = d\ is attained. Going 
away from that maximum, r\{a) has to decrease (if A / 0), since nontrivial characters 
oscillate. From the iVb = N c = 1 example (and for reasons which, if not clear already, will 
soon become evident), we expect a finite first derivative and hence a cusp singularity in r\(a) 
at a = 1. The task is to compute the slope of r\(a) on the a < 1 side of that cusp. 

To that end, we start from the formula 

2 lim ^-r x (a) = - lim ± [ X x ( Det"^ (l - X 2 ) dX . 

u(iVc) 

We would like to take the derivative inside the integral but, as it stands, are not allowed to 
do so. Since the obstacle arises from the noncompactness of the integration domain, the trick 
will be to compactify it. 

In the first step, we introduce the eigenvalues of X and the diagonalizing matrix as the new 
variables of integration. Thus, let x = {x\, . . . ,xn c ) € (iR)^ be the set of eigenvalues of the 
anti-Hermitian matrix X, and let t + C t = (iR)^ be the positive Weyl chamber given by 
\x\ < 1x2 < • • • < ixn c . If T = U(l)^ c is a maximal torus, say the diagonal matrices in U(iV c ), 
some dense open set in u(N c ) is diffeomorphic to t + x \J(N C )/T by the polar coordinate map 
tp : (x,gT) i-> gxg~ x = X. The Jacobian of this map is A(kc) 2 = \N c {N c -\) Y\ k<l {xk - xi) 2 , 
which means there exists some positive flat density dx on t + such that 

tf;*(dX) = dg T ■ A(ix) 2 dx , 

with dgT an invariant volume form on U(iV c )/T. 

The next step is to write the character xx as a sum over the integer weight lattice L\ of the 
representation A: 

Xx(e W )= e iEn * e *- 

L\ is determined in principle by division of the elementary antisymmetric torus functions in 
Weyl's formula Xx( eie ) = £,x+p{ el9 ) /Cpi^ 1 ®)- Inserting the sum into the integral we obtain 

-'^s-^)=^wm^j E n(i^) d-4)-* ao,)^. 

t+ {n}Gi A k=l 

For further treatment, the restriction of the domain of integration to the positive Weyl cham- 
ber t + C t is inconvenient. Because the integrand is invariant w.r.t. the symmetric group (the 
Weyl group of U(iV c )), we may actually lift the restriction and integrate over all of t = (iR)^ c . 
The new integral is JV C ! times the old one, so we divide by that factor. 
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Now is the point where we compactify. Observe that 



/(*) = - 1 
JK ' l + z 

as a function of z G C satisfies \f(z)\ < 1 for Rez > 0, and < 1 for Rez < 0. 

Therefore, given some fixed term {n} in the sum over weights, we may modify the integration 
domain at infinity without changing the value of the integral, as follows: for all k with > 
we close the integration contour for the variable xt £ iM around the right half-plane Rez > 0, 
whereas for all k with n/% < we close around the left half-plane (for nj. = we may close 
either way). Having closed the contours for all the variables, we pull them in from infinity, 
by holomorphicity. The integration domain, say Cr n i, is now compact. 

After compactification, we are permitted to differentiate under the integral sign and set e to 
zero, which gives 



d N c 
lim_—r x (a)=N c r 1 vol(U(N c )/T) ]T ]Tn fc J x k - xf)'^ A(ix) 2 dx . 

{n}£L x k=l q 



{n} 

To compute the remaining integral, we first make the closed integration contours identical for 
all the variables x±, . . . , xn c . This is possible to arrange because after setting e to zero, the 
point at infinity is no longer a singularity for the variables x\ with I ^ k (although for I = k it 
still is, because of the presence of the factor x^ in the integrand). Thus we deform from Ci n \ 
to some C a (j lk ^ j x ... x C CT ( nfc ) — C^ n ^. The subscript u{nk) reminds us that C CT ( nfc ) lies in 

the right or left half-plane — and encircles the pole of (1 — x^)~ Nc at x^ = ±1 — depending 
on whether nf. is positive or negative, respectively. 

Next we symmetrize the integrand, replacing x k by iV^ 1 ^2iXi, and we revert to integrating 
over a restricted domain (and drop the factor l/iV c !), by requiring the variables x±,... ,xjy c 
to be arranged in ascending order on C CT ( nfc ) according to the orientation of C a r nk y 

Now let M CT ( nfc ) be the adjoint orbit of U(iV c ) on CP •.. Then, using the polar coordinate 
map in reverse we arrive at 

d , , ^ A f Nr lr TrXdX 



{n}£L x k=l M 



Det^(l -X 2 ) 

(»fc) 



This integral is easy to compute. The integration domain M ff ( nfc ) is a closed orientable N%— 
manifold in the complex space g[(N c ,C), and the integrand, a holomorphic density, can be 
viewed as a differential form (of degree N 2 ) which is closed. Therefore, the value of the 
integral does not change if we contract the integration manifold M a ^ nk ^ to another one in 
the same homology class enclosing the singular point X = sgn(nfc) x 1. At that singularity, 
the scalar factor N c ~ 1 TrX takes the value sgn(nfc). We extract the factor with this value 
from the integral. Having done so, we expand the integration manifold to its original form, 
and then see that the remaining integral is unity by normalization of the Cauchy distribution 
Det _Arc (l - X 2 ) dX. Thus we get the result 

lim — r\(a) = d\Css\{\) , 
a— >i- da 
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where 

N c 

Cas 1 (A) = ^ 1 £ El n ^l ( 31 ) 

{n}eL x k=l 

is some sort of Casimir invariant, as will be explained shortly. 
An equivalent way of stating the result is 

= d x (l-(l-a)Cas 1 (\)+K(a)) {a < 1) , 

c (a) 

with a remainder term lZ(a) that vanishes faster than 1 — a as a — ► 1, from which we deduce 

l im ( 1 ^(l-M/iVA " _ c - M Cas l( A) 

iV^oc V d A Co(l - /x/JV)/ 

Given this formula, the rest of the calculation goes the same way as in Section 14.31 and we 
just write down the answer. The continuum limit is defined by putting a p = 1 — A p /a 2 and 
taking the plaquette areas A p to zero. Then, setting /i = J2 p A p /a 2 we have the local data 

r(W,/i) = ^d x exp(-/iCasi(A)) X \{U) , (32) 

A 

and hence the genus g partition function 

A 

as follows again by the reasoning of |2*5] . 

What is the corresponding continuum field theory, i.e. what's the Lagrangian that gives rise 
to the Hamiltonian Casi upon canonical quantization? To answer that question, we first 
need to sharpen our understanding of Casi. Recall from Lie theory that Casimir invariants 
are elements in the center of the universal enveloping algebra (of the Lie algebra at hand), 
which consists of polynomials in the generators. The invariant Casi certainly does not come 
from a finite-degree polynomial, so we are not entitled to call it a Casimir invariant in the 
strict sense. However, it is something quite similar. Fix some basis {ta} of u{N c ) and write 
X £ u(iV c ) as X = ^2 a Xata- Reinterpret ta as an abstract operator ta that acts in the 
linear space V\ of any representation D x by D^(ta)- Then, viewing the coefficients Xa as 
real-valued functions on u(iV c ), consider the formal expression 



a; 



M J Bet Nc (fi 2 — X 2 ) ' 

u(7Vc) 

By the replacement ta — ► D^{ta), this makes sense as a compact operator in every unitary 
representation space V\. Moreover, by the invariance of the Cauchy distribution w.r.t. con- 
jugation X — > UXU~ 1 , this operator commutes with all the U(A C ) generators, so it is a 
multiple of unity in every irrep by Schur's lemma. We claim 

Casi = lim . 

M^0+ jji 
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Indeed, since acts as a multiple of the identity on every irreducible V\, we might as well 
take the trace over V\ and divide by the dimension d\. What we then encounter is an integral 
of the character xx against the Cauchy distribution, and by a slight variant of the calculation 
done earlier for r\(a), we find that ^/i|y has a right-hand derivative at \i = and the value 
of this derivative is Casi(A). 

When acting in the Hilbert space of square-integrable class functions f(U) on U(iV c ), the op- 
erator ^ x a t a caus es translations f(U) \— > f(Ue x ) = f(e~ x U). It follows that the invariant 
operator acting in the same Hilbert space has the integral kernel 

(U'\K^\U) = J 5 (U- x U'e x ) Det-^ (fi 2 - X 2 ) fi N " dX . 

u(JV c ) 

By definition, T(U~ 1 U' , /i) in (j32j) is the kernel of the one-parameter semigroup generated by 
Casi. From what we have just shown, the same holds true of (U' iTJjj) asymptotically for 
small /U. Hence we infer the limit relation 

r([/-V,/i) = Jim (u'\ (K^f \u) , 

which allows us to construct a functional integral representation for T(JJ, fi) in the stan- 
dard fashion. The resulting field theory is given by the action functional anticipated in the 
introductory part of the current section. Indeed, expressing the tangent-space J-function 
5(U~ l U'e x ) by Fourier integration over a conjugate variable (p with values in u(N c ), we 
encounter the Fourier transform of the Cauchy distribution: 

J e iTr0X Det -7V c ^2 _ X 2) ^ ^ 
u(N c ) 

On substituting X by \xX this becomes precisely the type of integral we computed in the 
current subsection, albeit with the character \\ ((1 — + eX)) replaced by the bounded 

function e lTr<J ^. Proceeding in the same manner as before we find that the result in the small- 
fi limit is approximated by e - ^"^" 1 with || 4> \\i the linear potential in (|14l) . 

4.5 A different perspective: Howe duality 

In treating the two-dimensional theory, crucial use was made of the expansion of the plaquette 
distribution function in terms of U(A r c ) characters x\(U)'- 

|Det (1 - a U) \~ 2Nb = cx(a) X x(U) . 

A 

We now wish to communicate the intriguing fact that the expansion coefficients c\{a) them- 
selves are characters, of a noncompact Lie group dual to U(-/V c ) in the sense of R. Howe 
|15[ I14j . Explaining this duality will go a long way toward preparing and setting up the 
duality transformation carried out in the next section. 

The general framework for the kind of duality we are about to describe is what is called the 
bosonic Fock space in physics (and the Shale-Weil, or metaplectic, or oscillator representation 
in mathematics). To get started, let there be a single oscillator or boson mode with operators 
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b, satisfying the canonical commutation relations [b, b<] = 1. There is a vacuum |0), which is 
annihilated by b, and the bosonic Fock space V is the linear span of the state vectors (6^) n |0) 
with n G N U {0}. The space V comes with a natural Hermitian scalar product, with respect 
to which b^ is the adjoint of b. In the presence of a chemical potential > 0, the partition 
sum of a single boson mode is 

oo 

Tr v e-^ tfe = ^ e_M " = (l - e^) -1 . 

n=0 

Now we enlarge the formalism by adding a color degree of freedom i = 1, 2, . . . , N c to the boson 
operators b l ,b 1 ^. (Of course the canonical commutation relations still hold, the Hermitian 
scalar product extends in the obvious manner, and the vacuum |0) is annihilated by each of 
the b 1 .) Let E kl be the N c x N c matrix with entry 1 at the intersection of the /c th row with 
the I th column, and zeroes everywhere else. To a "Hamiltonian" X = ^2 kl X kl E kl G u(-/V c ), 
we associate the partition sum 

Trye-ME^V'eE^^ = Det (l - e^*) -1 . 

Now we extend the formalism further by adding flavor a = 1,2, .. . , iVj,, and distinguishing 

between bosons of opposite U(l) charges, ±. Then, if N hos = £\ a (b J f b J f + tff V a ) is the 
total boson number, we get 

Tr v e-^-e^^ 6 ^ , 

which becomes the distribution function we have been working with on identifying e~^ = a 
and e x = U G U(iV c ). 
The linear mapping 

x^y.E xkl ( fo + at& + - b - a * b - 

a=l k,l=l 

is an isomorphism of Lie algebras (i.e. it preserves the commutator), and it is easily seen to 
exponentiate to a homomorphism of Lie groups, U i— > Tjj. Thus, U G \J(N C ) acts on the 
bosonic Fock space V by Tjj. Using Tjj we can write the previous formula in the concise form 

|Det (1 - a U)\~ 2Nb = Try cA os T v . 

This identifies the plaquette statistical weight function of the boson induced gauge model as 
a trace, or character, in Fock space. 

In addition to U(iV c ) there is a second group that naturally acts on the bosonic Fock space 
V. This is the noncompact Lie group U(A r b,A^b). To describe its action, we decompose Lie 

algebra elements Y G u(-/Vb, -/V^) into 7V b x blocks as Y = d)' with anti-Hermitian 

A = Y,A ab E ab , D = Y,D ab E ab , and with B = J2B ab E ab , C = J2C ab E ab being adjoints of 
each other; and we assign to Y an operator Y = — Y^ on V by 



f = J2 (A ah VfV+ b + B a V/V? f - C ab b j l% b - D^btb^ 

a,b=l j=l 
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The mapping Y i— > Y is an isomorphism of u(Ab, N^) Lie algebras. By exponentiating it, we 
get a unitary representation of U(Ab, Ab) on V. 

Thus we have two groups, U(iV c ) and U(A b , Ab), acting on the bosonic Fock space. Because 
U(A C ) acts on color and U(Ab,Ab) on flavor, the two group actions commute. This turns 
out to be a maximal property: you cannot enlarge either one of the two groups without 
compromising it. Put differently, U(Ab, Ab) is the centralizer of U(A C ) inside the big group 
of symplectic (actually, metaplectic) transformations of V. R. Howe calls such a pair of Lie 
groups a dual pair, see [Hj for a pedagogical introduction to the subject. 

Because U(A C ) is compact, general theory guarantees that we can decompose V into irre- 
ducible representation spaces V\ for this group. All U(iV c ) representations of a given type 
A are collected into what is called an isotypic component for U(A C ) in V. Here comes the 
main point of the present discussion: as a consequence of the dual pair property, the prod- 
uct U(iVb,iVb) x U(A C ) acts irreducibly Jl] on every such isotypic component. Thus the 
decomposition of V takes the form of a multiplicity-free sum, 

V = Y j V x ®V x , 

A 

where V\ and V\ are irreducible representation spaces for U(Ab, Ab) and \J(N C ), respectively, 
and the correspondence V\ <-> V\ is one-to-one. 

We can now use this to decompose the trace of the product a Nbos Tjj over V. The operator 
a Nhos is trivial on the second factor of any isotypic component V\® V\ in V, while Tjj is trivial 
on the first factor. As a result, the trace over every isotypic component separates into two 
factors, one depending on a and the other on U : 

Tr^ A a^Tu = Tr^aA- x Tv Vx Tu = c x (a) X x(U) . 

The factors are of a similar kind: both are primitive characters; X\(U) for U(A C ), and c\(a) 
for U(7V b , iVb). What we have learned then is this: the coefficient c\(a) in the character 
expansion of the induced statistical weight function, 

|Det(l-al0r 2Afb =Tr v a^Tu = J2cx(a) X x{U) , 

A 

is itself a character; it is the value on a Nboa of the U(Ab, Ab) character associated to the 
XJ(N C ) representation A by the dual pair correspondence. 

Precisely speaking, iVbos does not represent a generator of the real form u(iVb,iVb), but a 
generator of the complexified Lie algebra g[(2iVb,C). (The generator lying in u(iVb, Ab) is 
iAbos-) Thus a Nbos is to be viewed as representing an element of GL(2Ab,C), and the c\(a) 
are obtained from U(Ab, N^) characters by analytically continuing to that element. 

Because the representation spaces V\ are infinite-dimensional, the characters c\(a) = Try a Nb 
all diverge at the unit element, a = 1. One of the results we found earlier is co(a) ~ 
(1 — a)~ 2NhNc+Nc near that point. This has a transparent interpretation from the present 
perspective: taking the trace of a Nbos over all states of V would give (1 — a)~ 2NbNc , since 
there are 2N^N C boson species. The projection onto U(A C ) singlets (the trivial representation, 
A = 0) amounts to imposing A<? constraints, which reduces the degree of the singularity to 
2A b A c -A c 2 . 
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Every irreducible representation A of U(iV c ) need not occur in the sum V = J2\ ® V\- 
Howe's statement is that the sum is multiplicity-free, i.e. the multiplicity is at most one, but 
it can also be zero. For example, for N c = 2 and N\, = 1 all the representations A = (Ai, A2) 
with Ai > A2 > or > Ai > A2 are missing. For > N c , however, all irreducible 
representations of U(iV c ) do occur, as is implied by the 5-function property 

5 Duality Transformation 

Having established the continuum limit of the 2d induced gauge model (J3J) (with TVjt, > iV c + 1 
boson species) to lie in the universality class of 2d Yang-Mills theory, we certainly expect 
a similar scenario to be true in higher dimension. The simple reason is that going up in 
dimension enhances the collectivity of the fields, and thereby works in favor of universality. 

The 2d Cauchy-type theory we found for iVb = N c (Section l4.4j) is a low-dimensional gimmick, 
and is highly unlikely to persist above two space-time dimensions. Thus in three dimensions 
and higher, we expect universal Yang-Mills physics in the continuum limit of the boson 
induced gauge model for all > N c . 

Motivated by this expectation, we will show in the current section how to pass to a dual 
version of the induced gauge model. The main tool we are going to use is a variant of the 
color-flavor transformation I17j , which is based on Howe duality and the notion of dual 
pairs we have just sketched. 

As an aside, we mention that the same goal of constructing a dual theory was pursued in 
0, by similar techniques. However, the transformation carried out there does not extend to 
the good range N\, > N c (in fact, only the case = 1 was addressed in that reference) and 
therefore fails to be of relevance for Yang-Mills theory. 

5.1 Abelian duality (review) 

By way of preparation and for future reference, we first review the standard construction of 
a dual theory in the abelian case. The nonabelian duality transform described later on will 
be seen to reduce to the standard one in the abelian limit. 

Consider the partition function of the induced U(l) gauge model (jlj) with = 1 (and no 
fermions, Nf = 0): 

Z{a) = ! [dU]Y[\l-a p U(dp)\- 2 . 
J p 

As before, the model is placed on a (i-dimensional cell complex A with boundary operator d, 
and all of the real parameters a p are chosen close to (but less than) unity. To dualize such 
a model in a concise manner, we need a few basic facts |3f)j from discrete exterior calculus, 
which are rapidly summarized in the next paragraph. 

Recall that the boundary operator d is a linear mapping from the vector space of /c-chains into 
the vector space of (k — l)-chains, and the boundary of a boundary is always zero: d o d = 0. 
A &-chain b is called closed if db = 0. The Poincare lemma (when applicable) states that 
every closed /c-chain b is a boundary: b = dc. Objects dual to chains are called cochains: a k— 
cochain uj (the discrete version of a fc-form) is a linear function that assigns to every fc-chain 
c a real number, say {c,uj). The coboundary operator d : fc-cochains — > (fc+ l)-cochains (the 
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discrete version of the exterior derivative on forms) is denned by demanding that the discrete 
analog of Stokes' theorem be valid: (c, da;) = (3c, a;). 

Turning to the case at hand, we write the U(l) elements on links as U(l) = e l9 ( l \ where the 
gauge field 6(1) takes values in the additive group R/27rZ. We regard the 9(1) 's as constituting 
a 1-cochain 9, and may then write U(\) = e 1 ^ 1 ) = e 1 ^'^ . In terms of 9, the U(l) field strengths 
on plaquettes are U(dp) = e i<ap -^ = e i<p ' de> . 

Now, inserting the Fourier expansion of the induced statistical weight function, 

|1 _ ae i<?|-2 = (1 _ a 2yl £ e -MlB(l/aO e W ( 

into the formula for the partition function, we get 

Z(a) = " 4)^ [ [d9}J2e- MU+i{nAe) , 

p n 

where the sum is over all 2-chains n = ^n p p (n p £ Z) on A. The notation || n \\ a means a 
weighted sum of absolute values: || n \\ a = Sp m (V a p) l n pl- 

The information contained in the partition function is somewhat limited, and to expose the 
full physics of the model we need to couple it to an external source, say by inserting a Wegner- 
Wilson electric current loop [3123]. Let C be any closed loop on A. Viewed as a 1-chain, C 
defines a linear function (C, 9) on 9; and this function is invariant under gauge transformations 
9 i ► 9 ~\~ df, by Stokes' theorem and the fact that C is closed. The Wegner-Wilson loop is 
W(C) = Ee 1 ^'^, where the expectation value is taken with respect to the statistical measure 
given by the partition function. 

Next we use Stokes' theorem to go from (n, d9) + (C, 9) to {dn + C,9). We then change the 
order of summation and integration and carry out the ^-integral, which results in the 1-chain 
dn + C being constrained to vanish: 



n:dn=—C 



It is in this dual form — as a sum over 2-chains n with boundary —C — that W(C) will 
be seen to emerge as a special case from the nonabelian duality transform described in the 
sequel; cf. Section EH1 

Let us quickly review how further analysis of the dual theory is carried out, and what is 
the physics to be expected. For brevity, we do this only for the special case of d = 3 + 1 
dimensions and lattices A with trivial homology groups. 

Choosing some surface S with boundary dS = C, one shifts n — > n — S in order for the sum to 
be over n with dn = 0. Since we are assuming trivial homology, the Poincare lemma applies 
and guarantees for every closed 2-chain n the existence of a 3-chain a such that n = da. Of 
course, a is not uniquely determined: if a has boundary n, then so does the gauge transform 
a + dip by any 4-chain (p. Thus, solving the constraint dn = by n = da, one arrives at a 
sum over gauge equivalence classes [a] = [a + dip] . For the Wegner-Wilson loop one obtains 

W(C) = e-" 5a - 5 H" / e" l|aa|U • (33) 
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At this point one usually transcribes the theory from A to the dual lattice A*. The transcrip- 
tion proceeds by a canonical isomorphism (known in the cohomological setting as Poincare 
duality) which turns fc-chains on A into (4 — /c)-cochains on A*. In particular, the 3-chains a 
on A become Z-valued 1-cochains a* on A*, and the boundary da becomes the coboundary 
da* . Passing to the dual lattice has the virtue of revealing the true meaning of the dual theory 
as a lattice gauge theory with gauge group Z. Other than that, the passage is really quite 
unnecessary, and one may as well continue to work on A as we do here. 

Finally, using Poisson summation one relaxes the values of a from Z to R, at the expense 
of inserting a factor e^ a ' m < and summing over all closed Z-valued 3-cochains m on A (they 
must be closed in order to be well-defined functions on gauge equivalence classes [a]). The 
idea behind this last step is that the M-valued gauge field a should become the dual of the 
"photon", and the closed 3-cochain m (which turns into a closed 1-chain on passing to the 
dual lattice A*) is to be interpreted as the world lines of magnetic monopoles. 

On a hypercubic lattice, and in the Villain approximation 

e -[|8o|la e -*£p (9a) p , 

the photon field a is now readily integrated out to produce an effective action for the magnetic 
monopole current m. The physics of the resulting model is readily understood jH] when t is 
sufficiently small: in that case the magnetic monopoles are bound in neutral clusters, leaving 
the system in a Coulomb phase with a massless photon. Rigorous mathematical control on 
this scenario has been achieved in |1U1 lllj . 

In the present model, that conclusion is less immediate. At very short scales of a few lattice 
units the photon cannot be free, as the action || da \\ a is not quadratic but linear (and in 
fact nonanalytic, by the use of the absolute value)! However, we expect that a sequence of 
suitable real space renormalization group transformations will cause flow toward the quadratic 
action. (The renormalization process of thinning the degrees of freedom by summing over 
the short-distance fluctuations of n, should have a similar effect as taking convolutions of the 
nonanalytic distribution n \-t e - "™"", and a generalized version of the central limit theorem 
should take effect.) We note that the method of Frohlich and Spencer which does not rely 
on integrating out the "photon" to produce an effective action for the magnetic monopoles, 
looks promising for a check on this picture. 

In summary, our induced U(l) gauge model is definitely interacting at short distances, but 
Coulomb behavior and a free photon are expected to emerge at large scales in four dimensions 
(provided that the a p 's are close enough to unity). 

5.2 Nonabelian transform: one— link integral 

Our starting point is the partition function (QJ on any d-dimensional cell complex A. Although 
our approach is general and can handle fermions as well as bosons, we will restrict our attention 
to the case of bosons only. Switching the order of integrations, we will first do the gauge field 
integral and afterwards the auxiliary boson field integral. 

For a fixed configuration of the complex boson fields <p, the Boltzmann weight of the field 
theory partition function Q is a product over links. Hence, integrating over the gauge field 
amounts to doing a set of independent one-link integrals. To write down and compute the 
one-link integrals, we need a good notation. 
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pen + (l) qelT(l) 

Figure 4: The set of plaquettes adjacent to a link 1, 11(1), consists of two subsets denned by 
comparing orientations, which either match [II + (1)] or don't match [II - (1)]. 

Let 11(1) denote the set of oriented plaquettes p that contain the link 1 in their chain of 
boundary links. In formulas: p G 11(1) if dp = ±1+ . . .; see Figure@J The set 11(1) consists of 
two subsets: the plaquettes whose orientation agrees with that of 1 (dp = +1 + . . .), and those 
where the orientations disagree (dp = — 1 + . . .). We write the decomposition into subsets as 
11(1) = n + (l) U II - (1). Note that the cardinalities of the two sets Il + and II - are the same, 
since every plaquette occurs twice, once with each orientation. 

Now we fix some link 1, write £7 = £7(1) and II = II ± (1) for short, and denote the two 
boundary sites of 1 by n\ and 712 (in other words, n\ is the site where 1 begins, and n<i is 
where 1 ends). Then the gauge field integral pertaining to link 1 is 




Here we are using the same short-hand notation as in (J3J): what the quadratic expressions 
really mean is 

^ P (ni)LVp(n 2 ) = £ ^^(n 2 )[/ ! J<(n 1 ) . 

i,j=l a=l 

In what follows we trade the integral over £7 for a sum over dual degrees of freedom. 
5.3 Quantization 

To proceed, we need some definitions. The main idea is to carry out a kind of "quantization" 
(not to be confused with quantization in the usual sense of passing to a quantum Hamiltonian 
formulation of the theory): for every color i E {1, . . . ,N C }, every flavor a £ {1, . . . , Ny,}, and 
every oriented plaquette p G II = II + U II - , we introduce a boson annihilation operator bp 0, 

and the adjoint creation operator bp a \ These operators obey the usual boson commutation 
relations: 

(with the 6's and 6^'s commuting amongst themselves), and they act in a Fock space V with 
vacuum |0) characterized by 

f/>|0) = . 
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Elements of the gauge group, U E U(iV c ), act on V by unitary operators Tjj with the properties 
Tu\0) = \0) and 

= (pen+), (35) 

3 

Tub^rl = Y,< uji (q £n ")- ( 36 ) 

3 

Thus the 6 q for q E II - transform according to the vector representation of U(-/V c ), and the 6 p 
for p E II + according to the covector representation (£7 = f/ -lT ). For the creation operators 
$ the situation is reversed. 

Some mathematical background to this setup was given in Section [4.51 There, we explained 
how Ty arises by exponentiating an isomorphism of Lie algebras. For present purposes, just 
knowing the existence of Tjj and its specified properties will be sufficient. 

Consider now any one of the factors in the integrand of the one-link integral (|34|) , for p E II + . 
Using the operator T\j (or rather the part, Tj/(p), acting on the bosons associated with p) we 
can express it as a matrix element between coherent states (index summations suppressed!): 

e ^ P (n 2 ) ^Vp(«i) _ ^Q| e ^p(™2)6 P Tu(p) e b pVp( ni ) |o^) 

Verification is immediate by the XJ(N C ) invariance of the vacuum (Tu\0) = |0)), relation ()35j) 
and (0|e^ fe e fet< ^|0) = e^. Similarly, for q E II~, we have 

e <?q(ni)f/~Vq(™2) _ e V>q(n 2 )C7^q(ni) _ ^ Q | (n 2 ) feq ^ g&q^q ("1 ) | 

Here we used (U~ l ) % i = U Jl and the fact that switching from p £ II + to q £ II - interchanges 
the vector and covector representations of U(iV c ). Thus, although the matrix U was replaced 
by its complex conjugate U, the operator whose matrix element we take is still T\j. 



5.4 Projection on gauge singlets 

We return to the task of integrating over U. The benefit of "quantization" , i.e. interpreting 
the integrand of the one-link integral as the matrix element of an operator in Fock space V, 
is that all dependence on the gauge field matrix U now resides in 

Tu = n t u(p) n ^(q) , 

pen+ qen- 

which is just the operator satisfying the relations (|331 l3"fi|) . Thus, doing the [/-integral has 
become very easy; given the last two formulas in the preceding subsection, we immediately 
express the one-link integral as a coherent state matrix element of Po = J V ^ N \ Ty dU : 

I exp ^ (p p (n 2 )Uip p (ni) + ^ £q(ni)J7~Vq( n 2) I dU = 
\J(N C ) \pen+ qen- / 

exp ( Y ^p( n ^) b p+ Y ^q( n 2) & q) p o ex P ( Y 6 P M n i)+ Y 6 q^q( n ^ 
pen+ qen- pen+ qen- 



o;. 

(37) 
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What is the meaning of PqI According to (j^SJEH), the operator Ty acts on Fock space V by 
U(iV c ) rotations. If \ip) is any state in V, integrating the rotated state Tjj\ip) against Haar 
measure dU kills that part of the state which transforms nontrivially under U(iV c ), and leaves 
only the U(iV c ) invariant part. Hence Pq is simply the operator that projects on the U(iV c ) 
invariant sector of V! 

The present formalism would be useless if that sector — let it be denoted by Vo — were some 
complicated and obscure subspace of V. Fortunately, that's not the case and Vo has a very 
transparent description, as follows. 

Consider the set of all boson pair creation operators of the form 

N c 
i=l 

where p G n + , q G and a, b G {1, . . . , iVb}. Because the b % ^ transform as a U(iV c ) 

vector and the as a XJ(N C ) covector, the contraction formed by summing over equal 
colors i = j = 1, . . . , 7V C is a \J(N C ) scalar. Thus the E^ are U(iV c ) invariant, and so is every 
state created by acting with an arbitrary polynomial in these operators on the vacuum: 

K% E plt--- E p n t\°)- (39) 

It is a beautiful and powerful theorem of classical invariant theory ^5] that the states 1)39(1 
span Vo; i.e. every \](N C ) invariant state in V can be reached by repeatedly acting on the 
vacuum |0) with the pair creation operators E^ and taking linear combinations. 

An outline of the reasoning is as follows. Let r = x cardll + , with cardll + being the 
cardinality of the set II + (which is the same as the cardinality of II - ). For a ci-dimensional 
hypercubic lattice the value of r is 2(d — 1). The antihermitian operators built from the 
-Epq and their adjoints, taken together with all their commutators, span a Lie algebra that 
generates a unitary representation (on V with the usual Fock space scalar product) of the 
noncompact group U(r, r). Because the E^ are U(A r c ) invariant, it is clear that the Fock 
space actions of the two Lie groups \J(N C ) and U(r, r) commute with each other; they in fact 
centralize each other and constitute a dual pair in the sense of R. Howe. 

Recall from Section that the representation of a dual pair \J(N C ) x U(r, r) is irreducible on 
every \](N C ) isotypic component of V. In particular, this implies that U(r, r) acts irreducibly 
on the U(A r c ) invariant sector Vo of V. This irreducibility is just what we are claiming: every 
state in Vo can be reached by acting with the generators of U(r, r) on some U(iV c ) invariant 
reference state. (And, of course, if we take that reference state to be the vacuum |0), it suffices 
to act with the pair creation operators E^.) 

Let us summarize where we were prior to the digression characterizing the subspace Vo: we 
had identified the integral Po = Ju^ N \ T\j dU as the projector onto Vo, and we had expressed 
the one-link integral (|34| as a matrix element between coherent states; schematically, (|37|) is 
of the form 

I(ip,<p) = (p(n2)|Jb|p(ni)) . 
In the following we will think about this matrix element as a trace: 

I (ip, (p) = Try(-PoAp) = Tr Vo A v , with A v = |y(ni))(^(ra 2 )| ■ 
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We mention in passing that the color-flavor transformation in its standard form, introduced 
first in a supersymmetric setting in |16j , and applied to the present context in [5] , would result 
at the present stage of development if it was possible to express the trace over Vo as an integral 
over U(r, r) coherent states. However, this can only be done if r < N c /2, and fails for the 
range of r values of physical interest. (The problem is a problem of convergence: the U(r, r) 
representation spaces Vo for N c < 2r cannot be realized by square-integrable holomorphic 
sections; cf. the last part of Section H3) We must therefore proceed in a different manner. 

5.5 Integration over <p, tp 

We have written the integral 1)34(1 over a single matrix U = U(l) as the trace of some operator 
A v (l) in the U(iV c ) invariant sector Vo(l) of an auxiliary Fock space V(l). In the final step 
taken now, we multiply all one-link integrals together and arrive at an expression for the 
partition function (^Q) of the form 

Z = Tr Vo A , 

where Vo = ®i^o(l) is the tensor product of the U(iV c ) invariant spaces associated with all 
the links, and the operator A is the result of integrating J^A^l) over the boson fields <p, (p. 
It remains to describe A, which is what we do next. 

Recall from Section |2] that the boson p p hops from site to site of the chain of boundary 
links of the oriented plaquette p, and carries a color index i £ {1, . . . , iV c } and a flavor index 
a £ {1, . . . , iVb}. (p p enters through the mass term mt, iP of (j3J) and, by the manipulations 
of the previous subsection, it appears as a parameter in the operators A v = \p(ni))(p(n2)\ 
(schematic notation). We now write mb jP = m p for short. 

While the boson fields were originally coupled by the link matrices U(l), they are now com- 
pletely uncoupled (the mass term does not couple them, and the operators A v under the 
trace don't either.) Thus we can do the boson field integration for each component pp a (n) 
separately. So let us concentrate on some pp a (n) and calculate the corresponding integral. 

We start the calculation by noting that pp a (n) and its complex conjugate pp a (n) both occur 
exactly once as a coherent state parameter, multiplying some boson operator b or b^ in the 
exponents on the right-hand side of (|37|) . (The reason is that, since n is a site visited by the 
boundary chain of p, there exist exactly two links in dp that are attached to n.) To decide 
which they multiply — b, or b^ — , let li and I2 be those two links in the boundary chain of 
p that have the site n as a boundary point. Recalling that the lattice links 1 come with an a 
priori orientation determined by A, we are led to distinguish between three cases. 

1. Let the site n be the end point of li and the starting point of I2, and let the orientations 
of li and I2 agree with that of dp; see Figure 031). Then, from (|37|) we encounter the 
following integral: 

Tp{n) = J ^ p (n)# p (n)e- m " l ^ (n)|2 e b p (l2) ^ (n) |0)(0|e^ (n)f 'p (ll) . 

c 

(The color and flavor indices play no important role here and have been omitted.) Ex- 
panding the last two exponentials and doing the Gaussian integral we obtain 

00 ^ 

T p (n) = ^m; fc -(6 P (l 2 )) fc |0)(0|(6 p (l 1 )) fe . 

k=0 
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(i) (ii) (iii) 

li n 




Figure 5: For the purpose of doing the boson field integration, one must distinguish between 
the three types of corner shown here. 

Since (6^) fc |0)/\/M is a normalized state, we see that what the operator T p (n) does is to 
simply transfer bosons from link li to I2, while weighting each boson with the inverse of 
the mass m p . Note that this transfer happens for each color and flavor independently A 
second subcase of the present type occurs when the orientations of h and I2 are opposite 
to that of dp. In that case the variables <p p (n) and (p p {n) exchange their roles, but the 
operator that results on doing the integral is still the same. 

2. Next, let the site n be the starting point of both links li and I2, and let the orientation 
of dp agree (disagree) with that of li (resp. I2); see Figure [Sfii). Then from (|3*7|) we 
encounter the integral 

C p (n) = y"^ p (n)#p(n)e- m ^^ (n)l % 6 p (ll) ^ (n)+fe p (l2) ^ {n) |0)(0| 

C 

00 ^ 

= Y, m P k vS b i^K(h)) k \o)(o\. 

k=0 

Clearly, C p (n) creates (an indefinite number of) boson pairs in normalized states, weighted 
by powers of the inverse mass. Again, this happens for each color and flavor indepen- 
dently. If the orientation of p is reversed, (p p (n) and (f p {n) again switch roles but the 
form of the operator C p (n) remains unchanged. 

3. Finally, let the site n be the end point of both links li and I2 (and, although it makes no 
difference, let the orientation of dp agree with that of I2); see Figure Etiii) . Then from 
(|3*7)) we have 




C 



00 

= ^m p " fe |0)(0|(6p(l 1 )fep(l 2 )) fc A!. 

fc=o 

The operator D p (n) now annihilates pairs of bosons (again in normalized states, and 
weighted by powers of the inverse mass; and if the orientation of p is reversed, the final 
result does not change). 



33 



5.6 Summary: dual theory 

In summary, doing the integral over the boson fields (p, (p in the Fock space formalism produces 
operators that either transfer bosons (T) from one link to a neighboring one, or create/destroy 
boson pairs (C/D) on adjacent links. In all these processes, the plaquette label is conserved, 
and so are the color and flavor quantum numbers. (In the latter two processes, pairs of 
bosons are always created/annihilated in identical color and flavor states.) We are going to 
need names for the sets of transfer, creation and destruction sites on dp; let these be t p , c p , 
and Op, respectively. 

A simplification now comes from the fact that taking the trace picks out the boson number 
conserving processes. Using this to extract a common boson mass weight factor, we conclude 
that the partition function of the theory in the dual representation is 

/ „ N h n c \ 

z=Tr Vo nv^nn^ > ( 4o ) 

\±p a=li=l / 

where N p is the operator counting the total number of bosons associated with p, and the 
operators Ap a have the following structure: 

4 a = n c p» n t p» n w • 

nSCp ngtp nGDp 

The left and right operators create resp. destroy pairs of bosons in identical states, while 
the middle ones simply transfer bosons. (The operator Ap a is illustrated for the case of a 
triangular plaquette in Figure 0) The trace runs over the linear space of states Vo spanned 
by all polynomials 

N c 

*5&(li)*S&0W • • • KltMO) where < q (l) = £ ^(1)6^(1) , 

i=l 

with pj G n+(lj) and q., G n~(lj). 

While the evaluation of this partition function remains a challenging task in general, there 
are two conclusions we can draw immediately. To arrive at the first one, we note that one 
may employ any complete set of linearly independent states to compute the trace (|4()|). If 
the projection on the U(iV c ) invariant sector Vo were absent, we could use a basis of states 
labeled by occupation numbers {rip a (l)}, but since the projection acts on the color degrees of 
freedom, this is not possible unless iV c = 1. Nevertheless, the occupation numbers summed 
over colors, 

N c 

£<(l) = n£(l), 

i=l 

are still good quantum numbers. Now reconsider the triangular plaquette with corners 
ni,n2,Ti3 in Figure H3 Because the operator Tp' a (n2) simply transfers bosons without chang- 
ing any of their quantum numbers, the contribution of a state in Vo to the trace (|40|) vanishes 
unless 
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pair annihilation 




pair creation 



m,n2 



transfer 



Figure 6: Illustration of the operator Ap a for the case of a triangular plaquette: at site 77,3, 
bosons are annihilated in pairs; at site ri2, they are transferred according to the flow of the 
links; at site n±, boson pairs are created. 



Similarly, because C p ' a (ni) and D l p a {n^) create and destroy pairs of bosons in identical single- 
boson states, the contribution of a state to (|4Uj) vanishes unless 

n -pQ"n 1 ,n2) = n p (^1,713) j n -pQ'n 1 ,riz) = n pV"n,2,n3) ■ 

Altogether, this leads to the conclusion that the occupation numbers (summed over colors) 
must be the same for each of the links in the boundary chain of p: 



n p ,a (l) = ?ip (independent of 1) 



(41) 



By identical reasoning, this conclusion holds true not just for the triangular plaquette p in 
Figure El but for any plaquette on A. 

The second conclusion we can draw immediately results from the fact that all states in Vo 
are created by the U(A r c ) invariant pair operators £pq(l), where p £ n + (l) and q G n~(l). If 
we write h p = ^2^=1 n p f° r the occupation numbers summed over flavors (as well as colors), 
this implies 

P en+(i) qen-(i) 



for every link 1. This equation is easy to interpret: it is equivalent to the constraint on the field 
strength of the Z-valued gauge field dual to the U(l) gauge field in U(A r c ) = SU(A r c ) x U(l). 
We will come back to this in Section [5.81 



5.7 Wegner— Wilson loop 

Up to now we have concentrated on the partition function of the induced gauge model. To 
access its full physics, gauge-invariant correlation functions such as the Wegner- Wilson loop 
must be computed. 

A gauge-invariant correlation function that contains the Wegner- Wilson loop and neatly fits 
into our formalism is constructed as follows. Let C be any closed oriented contour of total 
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length L on the d-dimensional cell complex A. Viewing C as a one-chain, we can write it as 
a sum of oriented links: C = \\ + I2 + . . . + 1l, where the end point of is the starting point 
of We then define U(C) as the ordered product 

U{C) = U{\ L )U{\ L - 1 )---U{h)U{\ l ) , 

and consider the expectation value (. . .) of the gauge-invariant function Det(l — acU(C)) 
with respect to the statistical measure given by the partition function of the lattice gauge 
theory. This generates the Wegner-Wilson loop W{C) in the fundamental vector representa- 
tion: 

(Det(l -a c U(C))) = 1 - a c W{C) + 0{a 2 c ) , W{C) = (Tr U{C)) . 

We now set up the calculation of this correlation function, as follows. Let d\ k = Tik+i — n^, 
i.e. n^+i is the end point and n k the starting point of the link in the one-chain of C. 
Introducing complex fermion fields ipCi associated with the sites visited by the contour C, 
we modify the primordial action Q by adding 

L 

Sc-,a c = ^2 ($c(n k )^ c (n k ) - aJ/ L ^c(rafc+i)C/(lfc)V ; c(«fc) 

k=l 

By the standard rules of fermionic integration we then have 

J Q—Sb—Sc;a C 



(T>et(l-a U(C))] 



To pass to the dual description, we proceed in much the same manner as before. There are 
a few minor changes in order to take into account the fermions, but these do not affect the 
general strategy: 

The one-link integral (|34|) is modified for all links 1 in the contour C by the additional presence 
of the fermion fields. To carry it out, we again "quantize", i.e. we introduce colorful fermion 
operators /^(l) and (1) which satisfy the canonical anticommutation relations and act 
on Fock space with vacuum |0). The U(iV c ) transformation behavior is still determined by 
relative orientation, now between 1 and C: 

T u(l) f c (l)Tl m = fcW ji (l) if C = +1 + . . . , 

i 

T u(1) f c (l)Tl m = ^./y.il)r"il) if C = -!+.... 

3 

Coherent states work formally the same way for fermions and bosons; we can still express 
^Utp as a ma t r i x element of Tjj between coherent states e^^O). 

Doing the integral over U(l) with Haar measure, we still get a Fock space projection operator 
-Po(l)) although the U(iV c ) invariant subspace it projects onto, is enlarged for every link 1 in 
C: there exist additional states, created by boson-fermion pair operators 



^(1) = E -^(1)^(1) (qen-(l)) if C = +! + ..., or (43) 

i=l 

F p,cW = E^^W (P en+ (!)) if C = -l+..., (44) 
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and by invariant pairs involving two fermions. The latter, however, will play no role ultimately, 
as we only want W(C) in the vector representation (corresponding to a single-quark source), 
so we refrain from writing them down. The total set of U(iV c ) invariant pair operators together 
with all their (super-) commutators form a certain Lie superalgebra (of $jl type), 0(1). 

We denote the enlarged U(iV c ) invariant subspace associated with a link 1 on C by V (l). 
Howe's theorem j£] on the multiplicity-free action of dual pairs in Fock space [the dual 
pair now being U(iV c ) with g(l)] still holds in the present generalized setting, so V (l) is still 
obtained by acting with all U(A" C ) invariant pair operators on the vacuum |0). 

Integration over the fermion sources works the same way as for bosons (Section I5.5J1 . 
The final result of doing the integration (for fixed color quantum number i) is an operator 
which is denoted by A l c and is defined as follows. A l c is a product over all sites n visited 
by the contour C. In the same way as was shown in detail for the bosons in Section 15.51 
the definition of each factor A l c (n) depends on the a priori orientation of the links in C that 
begin or end on n. For sites n of type (i) (Figure EJ, Aq transfers fermions; for sites of type 
(ii) it creates pairs of fermions; and for sites of type (iii) it annihilates pairs of fermions. 

All this discussion is summarized in the following final formula: 

/ . N c N h \ 

<Det(l-a c ^(C))> = Z- 1 Tr Vo;C a^n^n^^n^ > ( 45 ) 

\ i=l ±p a=l / 

where Nf counts the total number of fermions, and Vo ; c is the tensor product of the spaces 
Vq(1) (resp. Vo(l)) for 1 contained (resp. not contained) in C. The normalization Z~ l is given 
by SU|). 

5.8 Abelian limit: G = U(l) 

In this subsection we briefly comment on the special case N c = 1, adopting the minimal model 
Nb = N c = 1. In that model, both the color and flavor degrees of freedom are absent and 
drastic simplifications occur. 

We begin with the expression for the partition function ()4Uj) as a trace over Vo- The evaluation 
of the trace is simplified by the observation that the states of Vo for N c = 1 are in one-to-one 
correspondence with sets of boson occupation numbers. Thus for = 1 they are labeled by 
n p (l) G NU{0}. By Eq. (jlTj) . these occupation numbers are 1-independent: n p (l) = n p . (We 
put the tilda here to reserve the name n p for another set of integer variables to be introduced 
presently.) They also satisfy the constraint (|4*2|) . 

The best way to deal with the constraint is to switch variables. Recall that in the full set of 
occupation numbers, {n p }, every 2-cell p occurs twice: once with its proper orientation (+p), 
and once with its orientation reversed (— p). We now switch variables from n± p 6 NU {0} to 

n p = n p — h-p G Z , and l p = rip + n_ p = |n p |, |n p | + 2, |n p | + 4, . . . . 

If we view p i— ► n p as a 2-chain n on the <i-dimensional complex A, the constraint (|42|) simply 
says that n is closed: dn = 0. Indeed, the question whether p belongs to II + (1) or II - (1), and 
hence what sign is attributed to n p in (j42[) . is decided by comparing orientations. This sign 
information is concisely encoded in the boundary operator d : 2-chains — > 1-chains. 
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The variables l p , on the other hand, decouple from the problem; the sum over each of them 
is a geometric series: 

m p p = (1 - a p ) aj, pl , 

« P e|n p |+2(NU{0}) 

where a p = l/mp p as before. The power L p appears because N p counts the total number 
of bosons associated with p and the L p links about p all have the same boson occupation 
number [Eq. (|41j) ]. Thus the partition function (|4()|) takes the final form 

P n : dn=0 P 

which coincides with the expression that resulted from the standard abelian duality transform; 
see Section 

Turning to W(C), we observe that the same argument that gave constant boson occupation 
numbers on plaquettes, leads to constant fermion occupation numbers along the contour C. 
To extract W{C) from the generating function (|45|) . we must place exactly one fermion on 
each link 1 in C . States in the U(l) invariant subspace Vq(1) are created by the invariant pair 
operators (|4^|44|) . which implies that the creation of a fermion along C is always accompanied 
by the creation of exactly one boson. Therefore, along C the relation 1)42(1 is modified to 

e+ tip- £ ^q = 0, 

P en+(i) qen-(i) 

where e = +1 if the orientation of C agrees with that of 1, and e = — 1 otherwise. If we 
switch from the variables n p ,n_ p to n p ,Z p , this constraint becomes dn = —C. Solving the 
constraint by setting n = da — S with dS = C, we obtain W(C) in the form 

6 Discussion 

The distinctive feature of the class of lattice models for U(A r c ) gluodynamics introduced here, 
is that they are induced from a pre— theory with A^,/f species of local bosons and/or fermions. 
The statistical measure of these lattice gauge models is a product over elementary plaquettes, 
with each factor being a ratio of determinants. 

The boson induced models have a critical point at unit mass, which allows a continuum limit 
to be taken. Our careful study in d = 1 + 1 dimensions showed that, if > N c + 1, this 
continuum limit is U(iV c ) Yang-Mills theory, with a specific ratio of the U(l) and SU(iV c ) 
couplings. The ratio goes to unity for iVb — > oo. In contrast, the continuum limit for = N c 
is not Yang-Mills but an unusual theory, which exists as a (super-)renormalizable quantum 
theory because the Cauchy distribution on u(iV c ) converges under convolution to a stable 
family of distributions distinct from the heat kernel family. 

Going up in dimension enhances the collectivity of the gauge field (by increasing the number 
of transverse gluons) and thus works in favor of "universality". Therefore, whenever our 
model induces U(iV c ) Yang-Mills theory in two dimensions, we expect it to do so in higher 
dimensions, at least generically. 

Although we concentrated on the special case of the gauge group being U(A r c ), a very similar 
treatment is possible for all classical compact Lie groups. Proposition (|12jl. which is the 
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mathematical basis for the existence of a continuum limit, carries over to the normalized 
distributions Det(l - aU)~ Nh dU on Sp(2iV c ) and SO(2iV c ), with thresholds N h > 2N C + 1 
and ./Vb > 2N C — 1, respectively. 

In the second part of the paper we subjected the boson induced gauge model to a duality 
transformation (a variant of the color-flavor transformation), which demonstrably reduces 
to standard abelian duality for G = U(l). The partition function of the theory in the dual 
formulation is the trace of a color- and flavor-diagonal operator Y\ Ap a acting on a tensor 
product of modules Vo(l) generated by quadratic U(iV c ) invariants on lattice links I. Again, 
the dual formulation is not restricted to U(N C ) but exists for other classical compact gauge 
groups as well. (This is because all of the compact Lie groups U(iV c ), Sp(2A r c ) and 0(N C ) are 
placed in Howe duality with corresponding families of noncompact Lie groups. For SU(iV c ) 
and SO(iV c ) the duality transform is more complicated EE 120] due to the existence of 
nonquadratic invariants of "baryon" type .) 

We do not know at present how far one can push the analysis of the dual theory. We believe 
it to be possible to develop a combinatorial approach to handle the operator Y\ Ap a , but the 
details have not been worked out yet. A more refined understanding of the asymptotics of 
the modules Vo(l) at large boson number is also needed. Although they are generated by 
quadratic invariants, they are not freely generated in the range iVjt, > N c , as is evident from 
the result c (a) ~ (1 - a)~ 2NbNc+N - [see for the U(iV b , N h ) character c (a) = Tr Vo cA 

We do not understand at present exactly how the Howe duality transform connects with other 
recent proposals, such as the quantum gravity formulation of SU(2)^ = 3 Yang-Mills theory in 
|13j . Also, to bring our induced gauge models closer to continuum gauge theories of current 
interest, one may ask whether they can be extended to models with robust supersymmetry on 
the lattice |31j . Another natural question to ask is whether one can construct a coherent -state 
resolution of the identity operator on the modules Vo, and by adding an extra time direction 
to four-dimensional space-time pass to a five-dimensional continuum theory (of topological 
gravity?) via the coherent-state path integral method. 

Let us finish by commenting on just the issue of existence of a coherent-state resolution. 
Consider the Howe pair U(iV c ) x U(Ab, N^) acting in a Fock space V with 2N^N C species of 
bosons, as described in Section 14.51 We are interested in the U(iV c ) invariant sector Vo- In 
particular, we would like a coherent -state integral representation for the invariant Hermitian 
scalar product on Vo- For N c > 27V b this is a standard problem with the following standard 
solution. Let M be the Hermitian symmetric space formed by dividing U(iVb,iVb) by its 
maximal compact subgroup. M is modeled by complex x matrices Z with noncompact 

domain Z^Z < 1; and g = ( J € U(iVb, N^) acts on the points Z of M by 



The coherent -state expression for any ip € Vo is a holomorphic section tp(Z), transforming 
under g G U(A^b, -/Vb) as 



g-Z= (DZ + C)(BZ + A) 



(V N *(g)<p) (Z) = Det-^(L> - ZB) <p {g' 1 ■ Z) . 



The invariant Hermitian scalar product of two sections (pi and y?2 is 





(1 - Z*Z)dZdZ , 



zzt<i 
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where dZdZ is a flat density (suitably normalized), and invariance means 

(<p 1 \ i p 2 ) = (V N °(g)<p 1 \V N °(g)<p 2 ) for g G U(iV b , A b ) . 

All this makes perfect sense as long as iV c is big enough (or A b small enough). However, 
when A c is decreased below 2A b , the density Det Nc ~ 2Nh (l - Z^Z) dZdZ b ecomes singular at 
the boundary of the symmetric domain Z^Z < 1, and the coherent-state expression for the 
Hermitian scalar product ceases to exist in the form given. This does not mean that it ceases 
to exist altogether. Indeed, for the case A b = N c we can easily see how to fix the problem. The 
singularity at the boundary indicates that the proper measure to use is concentrated in the 
boundary. Now, the boundary of the Z-model for M always contains the unitary group U(Ab) 
(the set of solutions of Z*Z = 1) as a U(A b , A b ) orbit, i.e. if U is in U(A b ), then so is its image 
g^ 1 -U = (D — UB)~ l (UA — C). A straightforward computation shows that the Haar measure 
dU on U(A b ) transforms under the action of g~ x as g~ l * {dU) = |Det(D - UB)\~ 2Nh dU. For 
the case A c = A b = A (and only in that case) the factor multiplying dU is canceled by the 
multiplier in the transformation law for sections <p(Z), so that the boundary integral 

(tpM = [ MU)ip 2 (U)dU 

JV(N) 

is invariant, and (by uniqueness) is the coherent -state integral representation of the Hermitian 
scalar product on Vo- This implies that, while the sections <p{z) are sections on the 2A 2 - 
dimensional symmetric domain M, the complete information about them in the special case 
A b = A c = N at hand is already encoded in the values they take on approaching the A 2 - 
dimensional part U(A) C DM of the boundary dM. 

We do not understand the details of the analogous construction of an invariant boundary 
measure for A b 7^ A c < 2A b , although we know on general grounds that such a construction 
must exist. Having a detailed understanding of that construction would enable us to extend 
the color-flavor transformation to the whole domain A c < 2A b . 



A Appendix: Calculations for the One— Plaquette Model 
A.l Fermion induced model 

We compute the expectation value of the Wilson loop, W(ctf, 0), for the fermion induced 
model on a lattice consisting of a single plaquette (Section 14.111 . 

Doing the same steps as in the proof of statement (|12j) . Section [3.31 we obtain 

W(a f ,0) = N c (e Wl ) ai /(l) a! , 

where 

(F) a = / F(e ie \ . . .) J] H " ^ 6j \ 2Nf II \ eWk - eWl \ 2 d9 ^ ■ ■ ■ d9 ^ ■ 

Putting 2;^ = e iefe , we rewrite this as a multiple complex contour integral over the unit circle 
U(l) C C : 

{F)* = {-^ f F( Z1 , = aZj)Nt " z <) 2d *i ■ ■ ■ dz Nc . (46) 

U(l)' 



Nc J'=l "3 k<l 
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The integrand is holomorphic on C\{0} for each of the variables Zf., which suggests evaluating 
the integral by contracting all integration contours to zero. There is a pole of order N c + Nf 
at zero in each variable of the normalization integral (F = 1). In the numerator {z\) a of the 
expectation value, the order of the pole in the distinguished variable z\ is reduced by one. 
Thus, by contracting the contours and evaluating the residues at zero, we obtain 



W(a { ,0)=N c (N c + N { -l)x 
where / is the function 

f(zi,...,z Nc ;a 



xN c +N{-2xN c +N(-l 



d^- 1 f{z l ,...,z Nc ;a i ) 



if 



■ ■■d. 



N c +Nf-1 

Z N C 



f(zi, 



,ZN c ;oif) 



= 2JV C =0 



k<l 



For small values of N c and Nf the derivatives are easy to evaluate using computer algebra. 
The results obtained in this way for iV c = 3 and Nf = 1, . . . , 6 were shown in Figure 



A. 2 Boson induced model 

Passing from the fermion induced model to its bosonic analog amounts to doing an analytic 
continuation, from positive integers Nf to negative integers — Nf,. Given Eq. (|46|) of the 
previous subsection, this continuation yields W(0,ah) = N c (zi) ab / (1) ab with 

/ f<3 (47) 

U(l)JVc J ■>' 

The strategy again is to contract all integration contours to zero. As compared with the 
fermionic case, we now encounter iV^-fold poles at Zj = in addition to the poles at Zj = 
that occur when N^ is less than N c . We separately consider the three cases Nf, > N c , N\, = N c , 
and N h < N c . 



> N c : In this case the integrand has an At, -fold pole at Zj = for j = 1, . . . , N c and 
no poles at Zj = 0. Application of the residue theorem yields 



W{0,a h ) = N c 



with / being the function 



d£- 1 ---d%- 1 z 1 f{z lt ...,z Ne -,a b ) 



d^-'-'-d^fizu... ,ZN c ;a h ) 



Zl =...=z N =a b 



N c 



f(zi, . . .,z Nc -a) = J 



7 N b -N, 



T^Jlizk-Zif 



k<l 



Again, for small values of N c and Ab we have employed computer algebra to evaluate the 
derivatives. The results for N c = 3 and Ny, = 1, . . . , 6 were shown in Figure |2 
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AT b = N c : As before, the integrand has apparent poles (now of order N\, = N c ) at Zj = 
for j = 1, . . . , N c and no poles at Zj = 0. The previous formula for W(0, a^) is still valid, but 
now it is possible to go further and get an answer in closed form. For that purpose we note 
that the polynomial nfc</( z fc — zi) has degree N C (N C — 1) and can be viewed as the square 
of a Vandermonde determinant, which in turn is a sum over permutations: 

W(z k - Zl f = j sgn(7r)(z 1 -a h r^- 1 ---(z Nc -a h r^~ 1 

k<l \ireS Nc 

= Yl S M™') (zi - a b ) w(1)+7r ' (1) - 2 • • • (z Nc ~ a b )-(^)+-'(^)- 2 . 

7r,7r'G5jv c 

Inserting this double sum into the integral we see that almost all of its terms give rise to an 
integrand which is actually holomorphic at z = (and hence holomorphic everywhere in the 
unit disc) for at least one of the integration variables. A true singularity in every one of the 
variables occurs only for the terms with equal exponents, 

±(z l -a h ) N c- 1 ---(z Nc -a h ) N <- 1 , 

which arise from the permutations with 7r(j) + ir'(j) = N c + 1 (for all j). There exist N c \ of 
such terms in the double sum, and these are the only ones that give a nonzero contribution 
to the integral. Looking at (|47|) we see that the resulting poles are simple in each variable, so 

<l) ab = (2^^1(1-^)-^ , 

and 

W(0,a h ) = N c (zx) ah /{l) ab = N c a h , 
which is the exact result quoted in the main text. 



iVb < N c : In this case there exist poles inside the unit circle both at Zj = «b an d at Zj = 0. 
Hence we must investigate the analytic structure of the integrand when q variables are sent 
to z = «b and N c — q variables to z = 0. The poles at these locations compete with the zeroes 
from the numerator nfc</(- 2 fc — z i) 2 ■ By power counting one easily sees that the integrand at 
the nominal singularities is actually holomorphic in at least one of the variables unless q = N^. 
In other words, nonzero contributions to the integral come only from sending variables 
to z = «b ? and the remaining N c — variables to z = 0. Adapting the Vandermonde 
determinant expansion of nfc<z( z fc — z i) 2 to such a location, say z± = . . . = zj^ h = «b and 
ZN b +i = . . . = zn c = 0, we see that the corresponding Laurent series starts as follows: 

LLj=l\ z 3 °bj Zj j=1 j=N h +l 

where the terms omitted are regular in at least one variable. 

Since there exist (^ c ) ways to divide N c objects into two sets of resp. N c — objects, 
and (^) x A^b!(Ac — A'b)! = N c l, the normalization integral is 

(l) Qb = (2^iV c !(l-^)^ . 
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In the numerator (zi) ah a nonzero contribution results only when the distinguished variable 
z\ is sent to z = (as opposed to z = 0). This happens in a fraction N^/N c of all cases, 
and so the Wilson loop expectation value is reduced by this very factor: 

W(0, at) = N c { -^± = N c ^a h = N h a h . 
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